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The four-dimensional Dirac-Schrodinger equation satisfied by quark-antiquark bound states is 
derived from Quantum Chromodynamics. Different from the Bethe-Salpeter equation, the equa- 
tion derived is a kind of first-order differential equations of Schrodinger-type in the position space. 
Especially, the interaction kernel in the equation is given by two different closed expressions. One 
expression which contains only a few types of Green's functions is derived with the aid of the equa- 
tions of motion satisfied by some kinds of Green's functions. Another expression which is represented 
\ in terms of the quark, antiquark and gluon propagators and some kinds of proper vertices is derived 

. by means of the technique of irreducible decomposition of Green's functions. The kernel derived 

' not only can easily be calculated by the perturbation method, but also provides a suitable basis for 

nonperturbative investigations. Furthermore, it is shown that the four-dimensinal Dirac-Schrodinger 
equation and its kernel can directly be reduced to rigorous three-dimensional forms in the equal-time 
Lorentz frame and the Dirac-Schrodinger equation can be reduced to an equivalent Pauli-Schrodinger 
equation which is represented in the Pauli spinor space. To show the applicability of the closed ex- 
pressions derived and to demonstrate the equivalence between the two different expressions of the 
kernel, the t-channel and s-channel one gluon exchange kernels are chosen as an example to show how 
they are derived from the closed expressions. In addition, the connection of the Dirac-Schrodinger 
^ ' equation with the Bethe-Salpeter equation is discussed. 

1. Introduction 

o . 

ly-j It is the common recognition that the Bethe-Salpeter (B-S) equation which was proposed early in Refs.[l, 2] is a 
• rigorous formalism for relativistic bound states. The prominent features of the equation are: (1) The equation is 
derived from the quantum field theory and hence set up on the firm dynamical basis; (2) The interaction kernel in 
^ • the equation contains all the interactions taking place in the bound states and therefore the equation provides a 
I ' possibility of exactly solving the problem of relativistic bound states; (3) The equation is elegantly formulated in a 
manifestly Lorentz-covariant form in the Minkowski space which allows us to discuss the equation in any coordinate 
frame. However, there are tremendous difficulties in practical applications of the equation, particularly, for solving 
the nuclear force in the nuclear physics and the quark confinement in hadron physics. One of the difficulties arises 
from the fact that the kernel in the equation was not given a closed form in the past. The kernel usually is defined 
as a sum of B-S (two-particle) irreducible Feynman diagrams each of which can only be individually determined by 
a perturbative calculation. This definition is, certainly, not suitable to investigate the subjects such as the nuclear 
I force and the quark confinement which must necessarily be solved by a nonperturbative method. This is why, as said 
in Ref. [3]," The Bethe-Salpeter equation has not led to a real breakthrough in our understanding of the quark-quark 
force" . Opposite to the conventional concept as commented in Ref. [4] that " The kernel K can not be given in closed 
form expression" , we have derived a closed expression of the B-S kernel for quark-antiquark bound states in a recent 
publication [5]. The expression derived contains only a few types of Green's functions which not only are easily 
calculated by the perturbation method, but also suitable to be investigated by a certain nonperturbation approach. 
Another difficulty of solving the B-S equation was attributed to the four-dimensional nature of the equation because 
the relative time (or the relative energy) would lead to unphysical solutions. Therefore, many efforts were made in 
the past to recast the four-dimensional equation in three-dimensional ones in either approximate manners or exact 
versions such as the instantaneous approximation [6], the quasipotential approach [7-12] and the equal-time formalism 
[13-16]. 

As one knows, the four-dimensionally covariant B-S equation for a two-fermion system is ordinarily formulated 
in a second-order differential equation with respect to the space-time variables in the position space. This kind of 
equation has been shown to have unphysical solutions with the negative norm. It was pointed out in Ref. [17] that" 
The appearance of the negative-norm B-S amplitude is a quite common phenomenon in the B-S equation" . A similar 
phenomenon was encountered in the Klein-Gordon (K-G) equation [18,19] 

(□, -f m2)^(x) = (1.1) 
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which was originally viewed as the wave equation satisfied by the single free fermion states. It is well-known that 
the K-G equation, as a second-order differential equation, has a solution with negative probability. This is because 
the wave function of the equation is determined not only by its initial value V'(O), but also by the initial value of the 
time-differential^ \t=o which would possibly cause the solution to have negative probability [18, 19]. Nevertheless, 
the Dirac equation [18, 19] 

{id^ - m)V(a;) = (1.2) 

where dx=j^d^ has not the negative norm solution because it is a first-order differential equation. As widely 
recognized, the Dirac equation gives a correct description of the single free fermion states. 

Analogous to the case of single fermion, the relativistic states for a two-fermion system may also be formulated 
by a set of first-order differential equations just as the Hamilton equation in Mechanics and the Maxwell equation 
in Electrodynamics which are equivalent to the second-order differential equations, i.e. the Lagrange equation and 
the D'Alembert equation respectively. Motivated by this idea, it was proposed in the literature [20-28] that the 
quark- antiquark bound system may be described by two coupled Dirac equations which are constructed in accord 
with the Dirac's Hamiltonian constraint formalism [29] such that [23-25] 

[iSxi - mi - Vi{xi,X2)]tp{xi,X2) = (1.3) 



{id^^ -1712- V2{xi,X2)]tp{xi,X2) = (1.4) 

where ip{xi,X2) denotes the two-fermion wave function, mi and m2 are the masses of quark and antiquark, Vi and V2 
stand for the effective potentials which are determined by the requirement of satisfying the Lorentz-invariance, the 
charge conjugation symmetry and a certain constraint (or say, compatibility) conditions. With a constraint imposed 
on the relative time, the above equations will be reduced to a three-dimensional eigenvalue equation. 

As emphasized in the previous literature [23-25], Eqs. (1-3) and (1.4) are built up within the framework of 
relativistic quantum mechanics and the interaction potentials are given in a phenomenological way although they 
are inspired by the quantum field theory and, as demonstrated in Ref. [23], are linked with the corresponding B-S 
equation. Obviously, in order to understand the Dirac-typc equations for the two-fermion system more precisely, it is 
necessary to give such equations an extensive investigation and an exact formulation from the viewpoint of quantum 
field theory. This just is the purpose of this paper. In this paper, we limit ourself to discuss the quark and antiquark 
{qq) bound states. The results certainly suit to other two-fermion bound systems. First we derive two first-order 
differential equations for the quark-antiquark bound states from Quantum Chromodynamics (QCD) which describe 
the evolution of the bound state with the total (center of mass) time and the relative time respectivc;ly. These 
equations will be called Dirac-Schrodinger (D-S) equation because the Dirac equation is, in essence, the Schrodinger 
equation in the relativistic case which is identified with itself as the uniquely correct equation of describing the 
evolution of a quantum state with time in the quantum theory. Next, we concentrate our main attention on the 
interaction kernel appearing in the D-S equation. We are devoted to deriving a closed and explicit expression of 
the interaction kernel. The kernel will be derived by two different methods: one is to utilize equations of motion 
satisfied by the qq four-point Green's function and some other four-point Green's functions in which the gluon field 
is involved; another is to employ the technique of irreducible decomposition of the Green's functions involved in 
the D-S equation. The first method is similar to that proposed previously in Ref. [14]. The kernel derived by this 
method has a compact expression which contains only a few types of Green's functions. The kernel derived by the 
second method is expressed in terms of the quark, antiquark and gluon propagators and some kinds of three, four and 
five-line proper vertices and therefore exhibits a more specific structure of the kernel. Especially, the kernel derived 
can not only be easily calculated by the perturbation method, but also provides a suitable basis for nonperturbative 
investigations. The D-S equation and its interaction kernel mentioned above are Lorentz-covariant. We will show 
how this equation and its kernel are reduced to the exact three-dimensional forms given previously in Ref. [15] 
in the equal-time Lorentz frame. It is well-known that the D-S equation is represented in the Dirac spinor space. 
This equation actually is a coupled set of sixteen scalar equations. In practical applications, sometimes it is more 
convenient to reduce the D-S equation to the Pauli spinor space following the procedure proposed in Ref. [30] . By this 
procedure, we will obtain an equivalent Pauli-Schrodinger (P-S) equation represented in the Pauli spinor space from 
the D-S equation. In the P-S equation, the interaction Hamiltonian is explicitly given in a series expression which 
has an one-to-one correspondence with the perturbative expansion of the S-matrix. To illustrate the applicability 
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of the kernels derived and the equivalence between the aforementioned two differerent expressions of the kernel, we 
will show how the one-gluon exchange kernels in the D-S equation and the corresponding interaction Hamiltonian in 
the P-S equation can be derived from the closed expressions. Finally, we will discuss the relation between the D-S 
equation and the corresponding B-S equation. 

The remainder of this paper is arranged as follows. In section 2, we will first derive two Dirac-type equations 
satisfied by the qq four-point Green's function. Prom these equations, the D-S equations obeyed by the B-S amplitudes 
will be derived by making use of the Lehmann representation of the Green's function [31]. Then, we will show how 
the four-dimensional D-S equation is reduced to the three-dimensional one. In section 3, the first explicit expression 
of the interaction kernel in the D-S equation will be derived by virtue of the equations of motion satisfied by the 
Green's functions involved in the Dirac-like equations. And, it will be shown how the closed expression of the exact 
three-dimensional kernel can be written out from the four-dimensional one. In section 4, the second expression of the 
interaction kernel will be derived by means of the technique of irreducible decomposition of the Green's functions. 
In section 5, the D-S equation will be reduced to the corresponding P-S equation. Section 6 will be used to give 
a brief derivation and description of the one-gluon exchange kernels. The last section serve to discuss the relation 
between the D-S equation and the corresponding B-S equation and to make some remarks. In Appendix A, we will 
describe the derivation of the equations of motion satisfied by the Green's functions which are necessary to be used 
in the derivation of the D-S equation and its interaction kernel. In Appendix B, the irreducible decomposition of the 
relevant Green's functions will be performed for the purpose of deriving the second expression of the kernel. 



2. Derivation of the Dirac-Schrodinger equation 



The Dirac-Schrodinger (D-S) equation satisfied by the qq bound states may be derived from the corresponding 
equation for the qq four-point Green's function which is defined in the Heisenberg picture as follows [32] 



g{xiX2;yi,y2)afSpa = (0+ T{N[MxiW0{x2)]N[iPp{yi)ii"M]} 0") 



(2.1) 



where tp{x) and il)'^{x) are the quark and antiquark field operators respectively, tpi^) ^iid ip {x) are their correspond- 
ing Dirac conjugates [18] 



^f{x) = C4> {x),'iP (x) = -iP'^ix)C- 



(2.2) 



here C = 17^7" is the charge conjugation operator, | 0"*=) denote the physical vacuum states, T symbolizes the 
time-ordering product and N designates the normal product which is defined by 



iV[Va(a;i)V|(x2)] = T[Va(a;i)V|(a;2)] - (0+ |T[Va(a;i)V|(a;2)]| 0") 



(2.3) 



It is emphasized here that the above normal product can only be viewed as a definition in the Heisenberg picture. 
With the definition shown in Eq. (2.3), the Green's function in Eq. (2.1) may be represented as 



Q{xiX2;yi,y2)aPpa- = G{xiX2;yi,y2)al3pa + Sp{xi - X2)a0Sp{yi -y2)pa 

where 

G{xiX2;yi,y2)afjpa = (0+ T{tl)a{xi)i)l{x2)ii) p{yi)'ipl{y2)} 0") 
is the ordinary qq four-point Green's function [18], 

S*p{xi ~ X2)o.p = 1(0+ |r{V'a(xi)V'^(a;2)}| 0-) 

= Sf{xi - X2)af{C-'^)^i3 = S'}r{X2 - Xi)0xCxa 



and 



SUyi - y2)p. = i(0+ T{^,(?yi)^:(2/2)} 0-) 

= CaTSF{y2 - yi)rp = {C~^) psS'jriyi - y2)sa- 



(2.4) 
(2.5) 

(2.6) 
(2.7) 



3 



in which 



Sf{xi - X2)ay = (0+ \T{iP^{xi)i,^{x2)}\ 0") (2.8) 

and 

SUvi - 2/2)5. = (0+ \T{rs{yi)i^l{y2)}\ 0-) (2.9) 

are the ordinary quark and antiquark propagators respectively [18]. It is clear that the propagators defined in Eqs. 
(2.6) and (2.7) are nonzero only for the quark and the antiquark which arc of the same flavor. For the quark and 
antiquark of different flavors, the Green's function defined in Eq. (2.1) is reduced to the ordinary form shown in 
Eq. (2.5) since the second term on the right hand side (RHS) of Eq. (2.4) vanishes. In the case of the quark and 
antiquark of the same flavor, the normal product in Eq. (2.1) plays a role of excluding the contraction between the 
quark field and the antiquark one from the Green's function. Physically, this avoids the qq annihilation to break 
stability of a bound state. It would be pointed out that use of tp^{x) other than tp{x) to represent the antiquark field 
in this paper has an advantage that the antiquark field would behave as a quark one in the D-S equation so that the 
quark-antiquark equation formally is the same as the corresponding two-quark equation in the case that the quark 
and antiquark have different flavors. 

The equations of motion which describe the variation of the qq four-point Green's function G{xi^X2;yi,y2) with 
the coordinates Xi and X2 may easily be derived from the QCD generating functional as described in Appendix A. 
The results are 

{idx-, - mi)a-/G{xj^^X2;yi,y2)jl3pa = SapS^{xi - yi)S'^{x2 - y2)/3a 

+Cai3^'^{x-i - X2)S*p{yi - y2)pa - {V''^') ^^Gl^^xi | xi,X2;yi,y2)j0pa 



{idx2 - rn2)i3\G{xiX2;yi,y2)a\pa = SpaS'^ix2 - y2)SF{xi - y\)c 



+Ca/3^ (Xi -X2)Sp{yi -y2)pa- - (T )i3xGl{X2 I Xi,X2;yi,y2)aXpa 



in which 



-la A" 



(2.11) 



(2.12) 



where g is the coupling constant, = ^ and T = — A''*/2 are the quark and antiquark color matrices respectively. 



G»(.x. 
0+ 



Xl,X2;yi,y2)a/3pa 

T{Al (x. )^a {x^)r0 {X2)i^p (2/1 )X (2/2) } 



0- 



(2.13) 



with i = 1,2 are the new four-point Green's function including a gluon field in it and the propagators were defined 
in Eqs. (2.6)-(2.9). It would be noted here that the terms related to Sp{yi — 2/2) in Eqs. (2.10) and (2.11) are 
absent when the quark and the antiquark have different flavors. The equations of motion satisfied by the Green's 
function defined in Eq. (2.1) may be found by substituting the relation in Eq. (2.4) into Eqs. (2.10) and (2.11) and 
by making use of the following equations as mentioned in Appendix A 



[{idxi - niiJa-ySpixi - X2)^(3 = -Cap6^(,Xi - X2) - (r"^)a7A^*(xi | Xi, 2:2)^/3 



[{idx2 - m2)0\Sp{x-i - X2)a\ = -Ca0S^{xi - X2) - {T^'')i}\Al*{X2 I Xi,X2)a\ 



(2.14) 



where 



A^*(xi I xi,X2U = t(0+ |T{A^(a;i)Va(a;i)V|(a;2)}| 0") 

with i = 1,2 are a kind of quark-antiquark-gluon Green's function. The results are 

{idxi - mi)a-fGji3pa{xiX2;yi,y2) = Sc,pS^{xi - yi)S^{x2 - y2)f}a 
-{^°''^)a'iQp{xi I xi,a;2; 2/1, 2/2)7/3^(7 



(2.15) 



(2.16) 
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and 

{id^2 - m2)(j\QaXpa{xiX2;yi,y2) = S0aS^{x2 - y2)SF{xi - yi) 

-(r '')0\Qt{X2 I Xi,X2;yi,y2)a\pa 

where 



^bf- 1 - - ^ (2.17) 



Q^{Xi I Xi,X2;yi,y2)a0pa- = {O+\T[{N[A-ixi)Mx,)r0ix2)]N[i;^{yi)£{y2)]} 

= C^iXi I Xi,X2;yi,y2)al3pa + ^1,*{Xi I Xi,X2)al3S*Fiyi -y2)pa- 



(2.18) 



here i = 1,2. In the above, the normal products defined in the same way as that in Eq. (2.3). In comparison of Eqs. 
(2.16) and (2.17) with Eqs. (2.10) and (2.11), wc sec, the terms related to S*p{yi - ^2) in Eqs. (2.10) and (2.11) 
disappear in Eqs. (2.16) and (2.17). Therefore, the equations (2.16) and (2.17) formally are the same as given in the 
case that the quark and the antiquark are of different flavors. 

Multiplying Eqs. (2.16) and (2.17) with the matrices 7° and 72 respectively, we have 



[*Jr ~ h^^H^)]a-yGjl3p'7{xi,X2;yi,y2) 

= {7i)ap5'^{xi - yi)S'ir{x2 - y2)/3a + Q)^pp„{xuX2;yi,y2) 



0\ _ n,.\QC (ry. „,^\„ J_rM) ^^-n,. „,^\ (2.19) 



,o,„ .4,.__.,_^c„..,_...^ ... ...... (2-20) 



and 

~ h^'^\x2)]pxQaXpa{xi,x2;yi,y2) 
= {12)0't5'^{x2 - y2)SF{xi - yi)ap + Qa^p^{xi,X2;yi,y2) 

where 

/iW(x|) = -m ■ + maf (2.21) 

is the i-th free fermion Hamiltonian, 



SaPpai^i I Xl,X2;yi,y2) = {fli^)a-fG'^{xi I Xi,X2;yi,y2)'yf3pa 
Sa0pa(x-^ I Xi,X2;yi,y2) = {^2')l3\Qt{X2 \ Xi, X2; yi, y2)aXpa 



(2.22) 

yaffpAX'i I xi,X2;yi,y2) = {''iY )f3\y^{X2 \ xi,X2;yi,y2)aXpa 
here 



f^r = -gih^n, ^2 = -9lh2T\ (2.23) 

As will be proved in section 4, the Green's functions Q1^{xi \ xi,X2;yi-,y2) are B-S reducible, therefore, we can 
write 

Q^''\xi\xi,X2;yi,y2) = j (fzi(i^Z2K^''\xi,X2;zi,Z2)Q{zi,Z2;yi,y2) (2.24) 

where K^^\xi,X2; Zi, Z2) {i = 1,2) are just the interaction kernels. With the expression given in the above, Eqs. 
(2.19) and (2.20) can be represented as 

[*Jr - 'T'''^K^)]a'rGifipa{xi,X2;yi,y2) = {li)apS^{xi -yi)S''p{x2 -y2)i3a .2 25) 

+ ] d^Zid'^Z2K'-'^\xi,X2;Zi,Z2)al3\rQxrpa{zi,Z2;yi,y2) 

and 

~ h('^\x^)h\Q<x\pa{xi,x2;yuy2) = (72)/?<T^^(a;2 - 2/2)S'F(a;i -yi)ccp ,^ 26) 

+ j d^Zid'^Z2K^'^\xi,X2\ZuZ2)al3\rG\Tpa{zi,Z2;yi,y2) 

From the above two equations, we may obtain two equivalent equations: one describes the evolution of the Green's 
function Q{zi, Z2;yi,y2) with the center of mass time; another describes the evolution of the Green's function with 
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the relative time. The first equation is given by summing up the two equations in Eqs. (2.25) and (2.26). Introducing 
the cluster coordinates 

X = 771X1 + 772X2, X = Xi - X2] 

y = Viyi + my2, y = yi - y2, (2.27) 
„. = ,i = 1,2 

" mi+m2' ' 

the equation may be represented in the matrix notation as 

[i§-^-Ho(^,l>f)]g{X -Y,x,y) , . 

= S{X -Y,x,y) + Jd'^Zd'^zK{X - Z,x,z)g{Z -Y,z,y) ^ ' 

where 

S{X - Y,x,y) = 5\x^ - yiHS'p{x2 - 2/2) + 5\x2 - y2HSF{xi - yi) (2.29) 

Hopl, If) = h^^^lti) + 2) (2.30) 

is the total free Hamiltonian, 
2 

K{X - Z,x,z) =J2^^'H^ - Z,x,z) (2.31) 

i=l 

is the total interaction kernel and t = Xa is the center of mass time. In the above, the translation-invariance of the 
Green's function and the interaction kernel has been considered. 

The second equation mentioned above is given by subtracting the equation in Eq. (2.26) with weight r]i from the 
equation in Eq. (2.25) with weight 772 



[i^-Ho{l,^)]g{X -Y,x^) 

= S{X -Y,x,y) + J d'^Zd'^zK{X - Z, x, z)g{Z - Y, z, y) 



(2.32) 

J — 1 ,u^, yj T J u zju ^xv — ^ >^\^ — ^ i '^i y ) 

where 

S{X - Y, x, y) = r,26\xi - yih^,SUx2 - 2/2) - ViS\x2 - y2)llSF{xx - yi) (2.33) 

Ho{l, t) = mh^^\-^i) - Vih^^K'^2) (2.34) 
which is the relative Hamiltonian, 

K{X - Z,x,z) = mK'^^HX - Z,x,z) - tjiK'-^\X - Z,x,z) (2.35) 

which is the relative kernel and t = xq is the relative time. 

By virtue of the well-known Lehmann representation of the Green's function ^(^i, ^;2; yi, 2/2) [31], one may derived 
the equations satisfied by the B-S amplitude from the above equations. The Lehmann representation as shown below 
can easily be written out by the procedure of inserting the complete set of the qq bound states into the Green's 
function ^/(zi, Z2; 2/i, 2/2) denoted in Eq. (2.1), then considering the translation-invariance property of the Green's 
function and finally employing the integral representation of the step function. 



g{X -Y,x,y) = EAl d^Q„e-^Q"(^-^) 

n 



^ 2w„ t Q»-w„+ie Qil+Wn-ie -» 



where 
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XQ„(X,x) = e-Q"^XQ„(x) = (0+I^^(xi)r(x2)}|n) 
XQ„(y,y) = e'«"^XQ„(y) = {n |a^{^(2/i)/(2/2)}| 0") 

are the B-S amplitudes describing the bound state and w„ is the energy of the state | n). Upon substituting Eq. 
(2.36) into Eqs. (2.28) and (2.32), then taking the limit: limQ0_K^^((5^ — w„) and finally performing the integration: 

/ d'^Ye~^^^ , one can find that 

[i^ - Ho(t, -^)]XP.{X, x) = J d^Yd^yK{X - Y, x, y)xp,{Y, y) (2.38) 

and 

[i^-HoQt,^)]XP,{X,x) = J d^Yd^yKiX - Y,x,y)xp,{Y,y) (2.39) 

where the subscript <r in the B-S amplitude designates the other quantum numbers of a bound state. In the above 

derivation, the fact that the functions S{X — Y,x,y) and S{X — Y,x,y) have no the bound state poles has been 
considered. Eqs. (2.38) and (2.39) arc just the wanted D-S equations satisfied by the B-S amplitudes. Eq. (2.38) 
describes the evolution of the qq bound state with the center of mass time t, while, Eq. (2.39) describes the evolution 
of the qq bound state with the relative time r. Clearly, both of the equations are all the first-order differential 
equations whose solutions are only determined by the initial amplitudes at the origin of times. 

Considering the translation-invariance of the B-S amplitude and the kernels as shown in Eq. (2.37) and in the 
following 

K{X -Y,x,y) = J 0e-^Q(^-^)if (Q, x, y) 

K{X-Y,x,y) = J^e-'Q(^-^'>K{Q,x,y) ^' ' 

one can obtain from Eqs. (2.38) and (2.39) the equations satisfied by the amplitude which describes the internal 
motion of the qq bound system 

[E - Ho(P, ■^)]xp,{x) = j d'yK{P, x, y)xpM (2-41) 



Ho{P,^)]XP,{x) = J d''yK{P,x,y)xp,{y) (2.42) 
Furthermore, in view of the Fourier transformation 



XP.{x)^j^.e-^''''XpM 

{2wy (27r)i' 

Eqs. (2.41) and (2.42) will immediately be transformed into the momentum space 



K{P,x,y) = J ^^e-Wfe.;f(p,g,fc), 



(2.43) 



[E - Ho(P, r)]XpM = J ^K{P, q, k)xpM (2.44) 
- H^{P, -^)]XpM = j (04^(^' «' k)xp,{k) (2.45) 



where 



(2.46) 
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in which 

h^'^ (pt) =pt -01 + mi7°, i = l,2, ^2 47) 

P =Pi+P2, q = mPi - V1P2, k = r]2qi -r]iQ2 

here qi and pi are the initial state and the final state momenta for i-th single particle. Both of the above equations 
are identified themselves with the eigenvalue equations for the qq bound system. 

The D-S equations given in Eqs. (2.41) and (2.42) or (2.44) and (2.45) are Lorontz-covariant. They allow us to 
investigate the qq bound states in any coordinate frame. Based on the space-like property of a bound state, it is 
admissible to discuss the bound state in a special equal-time Lorentz frame. In this frame, the qq four-point Green's 
function becomes 

g{xi^xi:yt,jB;ti-t2) = {Q+\T{N[^ixt,h)^%xi,h)]N[lp{yt,t2)i^"{M,t2)]}\0-) (2.48) 
The equation in Eq. (2.28) is now reduced to 

^ Ho(MM)]G(xi ^x2-,yi,y2;ti - ^2) 
= S{ti- t2)S(xi^X2; |7i, 2/2) + / d^zid^Z2dtzK{xi^X2; zt, zi; h - t^)g{zi^Z2; yt, yi; tz - h) 
where 

Sixt,x^; yt, = S^xt - yth^S^p{x^ -y^) + ^^(xt - y^HSpixt - yt) (2.50) 
in which 



(2.49) 



Spixt -ri) = i{0+ \T[i^ixt,h)[if;iri,ti)] 
SUx2 - = (0+ T[ip%x^,t2)i^'{y^,t2)] 



0-) 
0-) 



(2.51) 



are the equal-time quark and antiquark propagators respectively which are actually independent of time due to the 
translation-invariancc property. By the Fourier transformations 

g{xtr2; ti, h - 12) = |fe^^;(*i-*=)e(5t,a5|; tu yt; e) 

if(xt,5^;I?,^;ti-t2)=/+^|fe*^(*^-*^);^(xt,S5;2?,2?;i?) ^' ' 

Eq. (2.49) will be represented as [15] 

[E - HQ{xtx2)\Q{xi^x2;yt,y2;E) 53) 
= 5(^,^2; yi,y2) + J 2:1 (^^2:2-?i'(^,S2;^,^;-E)^(^_^;^,^;£) ' ' 

This just is the three-dimensional equation satisfied by the qq four-point Green's function defined in the equal- time 
Lorentz frame. 

In the equal-time frame, the relative time of the qq system is zero. Therefore, the equation in Eq. (2.32) becomes 
meaningless. We are left with only the equation given in Eq. (2.53). The Lehmann representation of the Green's 
function Q{xt^X2;yt,y2',t) is still represented in Eq. (2.36) except that the four- dimensional relative coordinates x 
and y in the B-S amplitudes are now replaced by the three-dimensional ones af and "y*. Substituting such a Lehmann 
representation into Eq. (2.53), by the same procedure as stated in Eqs. (2.36)-(2.38), one may obtain a D-S equation 
represented in the three-dimensional coordinate space such that 

[E - Ho{l, lt)]xp.{t, ^) = y d'^Yd'^yKil - 1^, ^, t)XP,{^, t) (2-54) 

where 

XP,(^,^) = e^^-^XP.(^) (2.55) 
Apparently, in the momentum space, Eq. (2.54) becomes [15] 

[E - Ho(P, t)]XP.{t) = J -0^K(P, t, t; E)xP,{t) (2.56) 

This is precisely the three-dimensional D-S equation satisfied by the amplitude xp^ ( ^ ) which describes the internal 
motion of the qq bound system and may be used to solve the eigenvalue problem for the system. This equation is 
rigorous because the retardation effect is completely included in the kernel of the equation. 
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3. Derivation of the interaction kernel 



In this section, the interaction kernel in the D-S equation wiU be derived by making use of equations of motion which 
describe the variation of the Green's functions with coordinates yi and 2/2- The equations satisfied by the Green's 
function G(xi,a;2; j/i, 1/2) are derived in Appendix A and can directly be written out from Eqs. (A. 22) and (A. 31) by 
setting the source J to be zero. The result is 

G{xiX2; yi, y2)al3Ta{i^yi + mi)rp = -5ap5'^{xi - yi)S''p{X2 - y2)(3a (-3 i\ 

-Cpa5\yi-y2)S*p{xi-X2)al3 + G't{yi\xi,X2;yi,y2)af}ra{X°'^)rp 



G{xiX2;yi,y2)a0ps{idy2 +m2)5a = -S0aS^{x2 -y2)SF{xi -yi] 



-CpaS^iyi -y2)Sp{xi - X2)al3 + Gl{y2 I Xi,X2;yi,y2)al3p5{T^'')5a 



(3.2) 



where G'^{yi \ x\,X2\ yi, 2/2) {i= 1, 2) were defined in Eq. (2.13) with the replacement of Xi by y^. 

Substituting the relation in Eq. (2.4) into the above two equations and employing the following equations obeyed 
by the propagator Sp{yi — 2/2) whose derivation was mentioned in Appendix A 



S'^piVi - y2)Taiidy, + mi)rp = CpaS'^iyi - 2/2) + aI* {yi \ j/i, y2)T<T(r'"')rp 



SpiVi - y2)p6{idy^ + 1712)50- = Cpa5'^{yi - 2/2) +A^ (2/2 | y\-,'yi)pb^ )«, 



where 



S^vixi 1 2/1, y2) 



it is easy to find 



= -(0+ 
I 



T{Kl{xi)^MWM} 



Q{X\X2\ yi,y2)aPTa{i^yi + n^l)rp = -SapS'''ixi - yi)S'jr{X2 - y2)l3a- 
+Q^{yi I Xi,X2;yi,y2)af3Ta{^'''")Tp 



and 



G{xiX2:yi,y2)ai3ps{idy2 + 111.2)50 = -Si3aS'^{x2 -y2)SF{xi -yi) 



ap 



+Ql{y2 I Xi,X2:,yi,y2)ocl3p5{Y"'^)5a 

where 



g«(2/i I xuX2;yi,y2)a0pa = (0+ T{N[i,^{xr)rp{x2)]N[Al{yi)i,MijM]} 

= Gt^{yi I Xi,X2;yi,y2)a0pa + Sp{xi - X2)al3A''* {yi I yi,y2)pa 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



here i = 1,2. 

To derive the interaction kernel, we also need equations obeyed by the Green's functions G'j^ixi \ xi, 2:2; 2/1, 2/2)- 
According to the description given in Appendix A, the equations for the Green's functions G^^^ixi \ a;i, 0:2; 2/1, 2/2) can 
be derived by differentiating Eqs. (A. 22) and (A. 31) with respect to the source J'^^{xi) and then setting J = 0. The 
result is 



G°(Xi I XiX2]yi,y2)al3Ta{idy^ +TOi)tp = -5ap5'^{xi -2/l)A^"(a;i I X2,2/2)/3o 
-CpaS'^iyi -2/2)A°*(Xi I Xi,X2)al3 + Gf^{Xi,yi I Xi,X2:,yi,y2)cxl3Ta{^°'^)Tp 



and 



G^(a;i I xiX2;yi,y2)ai3ps{idy2 + 1712)50- = -S(3^S'''{x2 - 2/2)AJ^(xi | xi,2;i)c 
-CpoS^{yi -2/2)A°*(a;i | xi,X2)a0 + G'^l{xi,y2 \ xi,X2;yi,y2)a0p5{r '')5„ 



(3.8) 



(3.9) 
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where 



A«(xi I xi,yi)^p = i(0+ |T[A«(a;0V7(a;i)Vv(2/i)]| 0-), 
A^*(xi I yi,2/2)xa = i (0+ |T[A^(xi)^.(yi)V^(y2)] 0") 



and 



(0+ T{A°(a;i)A^(y,-)V'«(a;i)VM^2)Vp(2/i)V'.(y2)} 0") 



(3.10) 



(3.11) 



here i,j = 1,2. On inserting the relation in Eq. (2.18) into Eqs. (3.8) and (3.9) and utihzing the equations in Eq. 
(3.3), we are led to 



G'HiXi I XiX2;yi,y2)af}Tcrii dyi+ mijrp 

+SfAxi,yi I a;i,a;2;yi,y2)a/3r<7(r"'')rp 



-SapS'^ixi - yi)K'j^{Xi I X2,y2)f}a 



and 



G'ilixi I xiX2;yi,y2)af3p6iidy^ +m2)so- = -^i3a^'^ix2 -2/2)A^(a;i | xi,yi)ap 



+Qnt{Xi,y2 I Xl,X2;yi,y2)af3p5iT )5a- 



(3.12) 



(3.13) 



where 



S^i^iiVj I a;i,a;2;2/i,y2)a/3p<T 

=(0+|T{A^[A«(a;0Va(a;i)VM^2)]A^[A^(%)V^p(2/i)Vv(2/2)]} 
= GfAxi,yj I a;i,a;2;?/i,y2)a/3p<T + A^*(a;i I a;i,X2)a/3A**(yj | 2/1,^2) 



0- 



(3.14) 



Now, let us multiply Eqs. (3.5) and (3.6) respectively with and 72 from the right and sum up the both equations 
thus obtained. In this way, writing in the matrix from, we obtain the following equation 



Q{xi,x2]yi,y2)[i-Q^i + %{yi,y2)] 
2 

= -S{xi,X2;yi,y2) - E G^'Hvi I xi,X2;yi,y2) 



(3.15) 



where 
d 



d d 

+ 



d 



dt' dyl dyl ayo ' 



(3.16) 



Ho{yi,y2) = h^^\yi) + h^^\y2) 
here h^^\yi) and h^'^\y2) were represented in Eq. (2.21), 

Q^^\yi I X\,X2;yi,y2)<xl3pa = G^{yi I Xi,X2;yi,y2)a0Ta{^T)Tp 
G^'^\y2 I Xi,X2;yi,y2)a0pa = Q^{y2 I Xi,X2;yi,y2)a0ps{^T)Sa 

in which 

-giiliT = -91 ill T , ^ ~9l2l2T =--9l2l2 T 



(3.17) 



(3.18) 



(3.19) 



and 'S'(xi, a;2; t/i, ^2) was defined in Eq. (2.29). 

Similarly, on multiplying Eqs. (3.12) and (3.13) respectively with 7^ and 72 from the right and summing up the 
both equations thus obtained, one can get 
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= -T6,1^°' {xi,X2]yi,y2)a0pa- - Qp^" (xi, X2; yi, y2)al3pa 



where 



n^il^°'{xi,X2;yi,y2)oci3pa = ^'^{xi -yi)(7?)„pAJ;''(a;i | X2,y2)i3a (3 21) 

+5'^{X2-yi){l2)l3T^l{Xi\x-^,yi)ap 



and 

Q^il^°'{xi,X2;yi,y2)a(3pa = Gll{xi,yi I Xi,X2;yi,y2)<xl3Ta{^^i)Tp 



(3.22) 



+ Qfu{Xi,y2 I Xi,X2;yi,y2)af3psi^7)sa 

With the definitions given in Eq. (2.22) and in the foUowing 

■R^'^{xi,X2;yi,y2) = f2r^«"(a;i, 0:2; 2/1,^2) (3.23) 

Q^'\x,,X2;yi,y2) = nrQ«"(xi,a;2;yi,y2) (3.24) 

wc can write from Eq. (3.20) the equation satisfied by the function Q^^^p^ixi \ Xi,X2;yi,y2)- In the matrix form, it 
reads 

g^'\xi \ xi,X2;yi,y2)[i-g^ + ^o{yi,y2)] (3.25) 
= -1l^'-\xi,X2;yi,y2) - Q^'-\xi,X2;yi,y2) 

Up to the present, we are ready to derive the interaction kernel. Acting on the both sides of Eq. (2.24) with the 

operator i— + ^0(2/1,2/2) from the right and employing Eqs. (3.15) and (3.25), we have 

J d'^zid'^z2K^'^^{xi,x2; zi, z2)S{zi, z2-,yi,y2) = Tl^'^Hxi, X2;yi,y2) + Q^''\xi,X2;yi,y2) 

-J d^zid'^Z2K'^''\xi,X2\zi,Z2) Yj S^'\yj I zi,Z2;yi,y2) i^-^Q) 

In order to obtain the kernel, we may operate on the above equation with the inverse of S{xi,X2;yi,y2) and the 
kernel on the RHS of Eq. (3.26) may be eliminated by the following relation given by acting on Eq. (2.24) with the 
inverse of the Green's function g{xi, X2', 2/1,2/2) 

ifW(a;i,a;2;2/i,2/2) = j d'^uid'^U2g^''\xi \ xi,X2;ui,U2)g~'^{ui,U2;yi,y2) (3.27) 

With these operations, one may derive from Eq.(3.26) a closed expression of the kernel K^^\xi, X2', yi, y2) such that 

K^'\xi,X2;yi,y2) = J d'^zid'^z2{'fi^'\xi,X2;zi,z2) + Q'-''\xi,X2;zi,Z2) (328) 
-'D^''\xi,X2;zi,Z2)}S^^{zi,Z2;yi,y2) 

where 

V^'\xi,X2;zi,Z2) 

= / n d'^Ukd*Vkg^''> {xi I xi,X2;ui,U2)g~^iui,U2;vi,V2) J2 S^^Kvj I 'i'l, 1^2; ^1, -^2) (3.29) 
fe=i j=i 

In the above derivation, existence of the inverses Q~^{ui,U2;zi,Z2) and S~^{zi, Z2;yi,y2) has been assumed. The 
rationality of the assumption will be illustrated later. 
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Clearly, the total interaction kernel appearing in Eq. (2.38) is given by the sum 

2 

K{xi,X2;yi,y2) = J2 K^'^'>{xi,x2;yi,y2) = / d*zid*z2{TZ{xi,X2; zi, Z2) 
+Q{xi,X2; zi, Z2) - 'D{xi,X2; zi, Z2)}S~^{zi,Z2; yi, 2/2) 



(3.30) 



where 

2 

Ji{xi,X2;zi,Z2) =^'R^'\xi,X2;zi,Z2) (3.31) 



i=l 



2 2 

Q{xuX2;zi,Z2) = 5^QW(a;i,a;2;^i,02) = ^ n'^^'G^^Kx,, z, \ xi,X2;zi,Z2)n^" (3.32) 

i=l ij=^ 

and 

T>{xi,X2;zi,Z2) 

= / n d'^ukd'^vk G^''\xi \xi,X2;ui,U2)G~^{ui,U2;vi,V2)G''^\yj \vi,V2;zi,Z2) i^-^^) 

k—l — 1 

We would like here to discuss the role played by the last term in Eq. (3.30). In view of the relation in Eq. (2.24) 
and the following relation 

G''^\yj \ xi,X2;yi,y2) = j d'^zid'^Z2G{xi,X2;zi,Z2)K''^\zi,Z2;yi,y2) (3.34) 
which also follows from the B-S reducibility of the Green's function C/^-'^ {yj \ xi,X2]yi,y2) and considering 

j d'^zid'^Z2G{xi,X2;zi,Z2)G~^{zi,Z2;yi,y2) = S^{xi -yi)S^{x2 -1/2) (3.35) 
the function 'D{xi,X2; zi, Z2) in Eq. (3.33) may be represented as 

'D{xi,X2;zi,Z2) = I '^d'^Ukd^VkK{xi,X2;ui,U2)G{ui,U2;vi,V2)K{vi,V2;zi,Z2) (3.36) 
■' fe=i 

which manifests itself the typical structure of the B-S reducibility of the kernel. Therefore, the last term in Eq. 
(3.30) just plays the role of cancelling the B-S reducible part of the remaining terms in Eq. (3.30). As a result of the 
cancellation, the interaction kernel given in Eq. (3.30) is really B-S irreducible, consistent with the ordinary concept 
for the kernel in a two-body relativistic equation. Inserting Eq. (3.36) into Eq. (3.30), writing in the operator form, 
we have 

KS = H+ Q-KGK (3.37) 

This can be regarded as the integral equation satisfied by the kernel K. 

Analogously, In accord with the definition in Eq. (2.35) and the expression in Eq. (3.28), one may write out an 
explicit expression of the kernel occurring in Eq. (2.39) 

K{xi,X2\yi,y2)_ = mI^''^\xi,X2-,y\-,y2) - mK^'^\xi,X2-,yi,y2) 

=J_d'^zid^Z2{Tl{xi,X2]Zi,Z2) + Q{xi,X2;zi,Z2) (3.38) 
-'D{xi,X2; zi,Z2)}S~'^{zi,Z2; yi, 1/2) 

in which 
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(3.39) 



where A stands for TZ, Q or V. 

Now, wc arc in a position to write out the interaction kernel appearing in the three-dimensional D-S equation. As 
demonstrated in Ref. [15], the kernel in Eq. (2.54) can be derived by the same procedure as for the four-dimensional 
kernel. The expression of the three-dimensional kernel shown below formally is the same as the four-dimensional one 
described in Eq. (3.30) except that it is now represented in the three-dimensional space. 



K{x{ x2]yiiy2]E) = J (fzi(fz2{'R{xi^xi\zi,zi) 

+ 2(^,^2; Zi, Z2\ E) - V{xi X2; Zi, Z2\ E)]S~'^[z{, Z2\ It 1,1! 2) 



(3.40) 



where TZ{xi^X2;zi,Z2),Q{xt^X2;zt,Z2',E) and T>{xt^X2;zt,Z2',E) can be written out from Eqs. (3.31)-(3.33) as 
follows: 



(3.41) 



in which 

7lW«(xt,^; zt, n) = 6^{xt - zthi^l"{xt | ^) + S^i^ - ^)72 A^(St | xt, zt) 



(3.42) 



here A^(x| | xt,yt) and A^(xi \ X2, 2/2) ^i''^ defined as in Eq. (3.10) except that the time variables in all the field 
operators are now taken to be the same and therefore they are time-independent due to the translation-invariance 
property of the Green's functions, 



Q{xt,x^;zt,z^;E) = ^ | Sf,^; ^, ^; E)^/ 

in which Q^ti^t, zj \ xt^X2; "?2; E) is the Fourier transform of the Green's function defined by 
^^^(^) ^ I xt^X2; ^1, ^2; h - t2) 

= (0+ T{7V[A°(x|,ii)^(St,ti)V'^(55,ti)]A^[A^(^,i2)?^(^,t2)V^'(^,t2)]} 0") 



(3.43) 



(3.44) 



and 



2 2 

V{xt,xi;zt,z^;E) = / n d^Ukd^Vk E ^T^l^^'^i^i I xt,xi;ut,ui; E) 
fe=i ',j=i 



xG ^iui^U2;lfi,lfr,E)gl'^''{zj \ lfi,lf2;zt,Z2;E)Q: 



(3.45) 



in which gj^^°'{xi \ Xi^X2;ui,U2; E) and gi''^''{'Zj \ 'v'i,lf2;'zi,'Z2; E) are the Fourier transforms of the following 
Green's functions 



gj^'^'^ixi I xt^X2;ui,U2;ti -^2) 
= (0+ T{Ar[A«(S|,ii)V(^,ti)V^(s5,ti)]iV[V^(«i,t2)V^'(n2,i2)]} 0") 



(3.46) 



and 

Qi'^^izj I iti,ii'2;zi,z2;ti - t2) 

= (0+ |T{iV[V'(iJ?,ti)V^(t}5,ti)]iV[At(^,i2)?(^,t2)?'(l5,f2)]}|0-) 
and 5~^(ii^i2; yt, y^) is the inverse of the function in Eq. (2.50). 



(3.47) 
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4. Another derivation of the interaction kernel 



The aim of this section is to give a different expression of the interaction kernel in Eqs. (2.38) and (2.39) which 
will be derived by means of the irreducible decomposition of the Green's functions Gf^ixi \ xi, X2',yi,y2) denoted in 
Eq. (2.18). First, we start from the relation between the full qq four-point Green's function G{x\X2;yi,y2) and its 
connected one Gc{xiX2;yi,y2) which is derived in the beginning of Appendix B [18, 33] 

G{xiX2;yi,y2)_= Gc{xiX2;yi,y2) + Sf{xi -yi)S'jr{x2 -^2) 

-sUxi-x2)s*Ayi-y2) ^ ' 

where all the fermion propagators were defined in Eqs. (2.6)-(2.9). Substituting Eq. (4.1) into Eq. (2.4), we get 

Q{xiX2;vi,y2) = Gc{xiX2-,yi,y2) + Sf{xi - yi)Sp{x2 - ^2) (4.2) 

where the last term is the unconnected part of the function Q{xxX2; yi, 2/2)- 

Prom Eq. (B.5) given in Appendix B, we obtain by setting the source J = that [18, 33] 

G^(a;i I xiX2;yi,y2) = Gl^{xi \ xiX2;yi,y2) +J^^{xi \ xi;yi)S''p{xi - yi)_^^ ,^ 
+Sf{xi - yi)A^^°'{xi I 3:2; ^2) - A°*(a;i | xi,X2)S*p{yi - ^2) - S^^ixi - X2)A.l*{xi \ yi,t/2) 

where i = 1,2, G'^^{xi \ x\X2'-,y\,y2) is the connected part of the Green's function G'^{xi \ x\X2'-,yi,y2) and 
A^*(xi I xi^X2)^ A^(a:i | xi;yi), A^"(xi | 2:2; 2/2) and A° (xj | 2/1,2/2) are the three-point Green's functions which are 
given in Eqs. (2.15) and (3.10). On inserting the decomposition in Eq. (4.3) into Eq. (2.18), we see that the last 
unconnected term in Eq. (2.18) is cancelled out. Thus, we have 

g'^{xi I xi,X2;yi,y2) = G'^^{xi \ xi X2; j/i, 2/2) + A^(xj | xi]yi)S'p{x2 -2/2) . . 

+Sp{xi - 2/i)A^"(a;i | X2;y2) - S*p{xi - X2)1^*{xi \ 2/1,2/2) 

where i = 1,2. Substituting the above expression in Eq. (2.22), we will obtain the decomposition of the function 

^W(xi I xi,a;2;2/i,2/2)- 

In the following, we are devoted to analyzing the terms on the RHS of Eq. (4.4) through the one-particle-irreducible 
decompositions of the connected Green's functions included in those terms. The decompositions have been carried 
out in Appendix B. 



A. The t-channel one-gluon exchange kernel 

First we focus our attention on the second and third terms in Eq. (4.4). According to the decomposition in Eq. 

(B.15), the three-point gluon-quark Green's functions A^(xj | Xj;yk) which is fully connected can be represented in 

the form 

A^(a;i I Xj;yk) = j d^ziE^^ixi \ Xj;zi)Sp{zi - yu) (4.5) 
where 

E;(x, I X,-Zi) = j d^Uid''u2Ali{xi - Ui)Sp{Xj - U2)r*"(wi I U2,Zi) (4.6) 

in which 

A-iix, - uj) = \ (0+ |T[A«(xOAt(«,)]| 0-) (4.7) 

is the exact gluon propagator and Y^" [u\ \ U2, zi) is the gluon-quark three-line proper vertex as defined in Eq. (B.17). 
The one- particle- irreducible decompositions of the three-point gluon- antiquark Green's functions A.^{xi \ Xj;yk) can 
be obtained from Eqs. (4.5) and (4.6) by the charge conjugation transformation. The result is 
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(4.8) 



where 



T,l^Xi I Xj;z2) = J d\id\2Af,{xi - ui)SUxj - «2)r^''(«i | U2, Z2) (4.9) 

in which T'^'^{ui \ U2,Z2) is the gluon-antiquark three-hne proper vertex as defined in Eq. (B.18). 

When we set i = j in Eqs. (4.5) and (4.8), it is found that the S^(a;i | xi;zi) in Eq. (4.6) and the ^^"(0:2 | X2; Z2) 
in Eq. (4.9) are respectively related to the quark and antiquark self-energies in such a way 

OrS^(a;i I xr,zr) = -7?S(a;i, ^1) = S(a;i, ^1), 

07S^«(a;2 I X2; Z2) = -1^2^%X2, Z2) = S^(a;2, ^2) ^ ' ' 

Thus, from Eqs. (4.5), (4.8) and (4.10), we have 

n^klix^ I xi; yi)S%{x2 - V2) + ir^''h.f{x2 \ X2; z)Sf{x, - y{) 

= Jd'^zid'^Z2'S{xi,X2;zi,Z2)SF{zi-yi)S^{z2-y2) ^' ' 

where 

Y,{X\,X2\ Zi,Z2) = Tj{xi,Zi)5'^{x2 - Z2) + YF-(x2tZ2)^'^[xx - Zi) (4-12) 

which is the total self-energy of the gg system. According to the definitions given in Eqs. (4.4), (2.22)-(2.24), (2.31) 
and (4.2), we see, S(a;i, 0:2; z\,Z2) as a self energy term to appear in the interaction kernel. 
In the case of i 7^ j, from Eqs. (4.5) and (4.8), it can be written 

A^(X2 I X,-y,)S%{X2 - y2) + '^T^l^^l I ^2; y2)SF{xy - yy) 
= J d'^Zid'^Z2Kt{xi,X2;Zi,Z2)SF{Zi - yi)S1r{Z2 - y2) ' 

where 

Kt{xuX2;Zy,Z2) = ^T^l{X2\xv,z{)5\x2-Z2) + nT^l%Xy \X2;Z2)5\xi-Zy) (4.14) 

Based on Eqs. (4.4), (2.22)-(2.24), (2.31) and (4.2), it is clear that the Kt{xi,X2;zi,Z2) acts as a part of the 
interaction kernel to appear in the D-S equation. As will be seen in section 6, this part is precisely the kernel of 
t-channel one-gluon exchange interaction.. 



B. The s-channel one-gluon exchange kernel 

Next, wc turn to the last term in Eq. (4.4). The one-particle irreducible decomposition of the three-point Green's 
function in this term may also be found from Eqs. (4.5) and (4.6) by the charge conjugation transformation. The 
result is 

A;*(xi I 2/1,2/2) = - j d''z,d''z2d^zA;i{xi-z)t''*{z I z,,Z2)SF{zi-y,)SUz2-y2) (4.15) 

where T^'^*{z \ Zi,Z2) is the gluon-quark-antiquark proper vertex defined in Eq. (B.19). With the decomposition 
given above, according to Eqs. (2.22)-(2.24) and (2.31), the contribution of the last term in Eq. (4.4) to the kernel 

K{xi,X2\yi, 2/2) can be found from the sum 

-^0,"^S^(a;i -a;2)A^*(a;i | j/1,2/2) = / d^zid''z2Ks{xi,X2;zuZ2)SF{zi - yi)S''F{z2 - y2) (4.16) 
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where 



K,{xl,X2■,z^,Z2) = S*Axi-X2) / ciS^f^rA;t(^. " ^W'^^z \ z^.z^) (4.17) 

j=i 

is just the s-channel one-gluon exchange kernel occurring in the D-S equation which will be discussed in section 6. It 

is noted here that the s-channcl onc-ghion exchange describes the annihilation and creation process which takes 
place between the quark (antiquark) in the initial state and the antiquark (quark) in the final state as indicated 
by the ghion propagator in Eq. (4.17) (not between the quark and the antiquark both of which are simultaneously 
related to the initial state or the final state B-S amplitude for a bound state). 



C. The kernel from the Green's function G^^^ixi \ xi,X2;yi,y2) 



Now let us concentrate our attention on the irreducible decomposition of the first term in Eq. (4.4). As stated in Ap- 
pendix B, this decomposition may be derived from the functional differential of the Green's function Gc{xiX2',yi,y2) 
with respect to the source J°^(a;j) by using the one-particle irreducible decomposition of the function Gc(xi^X2; yi, 2/2)- 
The latter decomposition whose derivation is sketched in Appendix B is well-known [18, 33] and can be represented 
in the form 

2 

Gc{xiX2;yi,y2) = I H d^Uid^ViSpixi - ui)S^{x2 - U2) 
xT{ui,U2;vi,V2)Sf{vi - yi)S^{v2 - 2/2) 
where 

3 

T{ui,U2;vi,V2) = ^I'i{ui,U2;vi,V2) (4.19) 
in which 

ri{u„U2;V„V2) = - I d^Z,d^Z2r'"'{z, I U„V,)Dl'i,{z, - Z2)r'^'''{z2 I U2,V2) (4.20) 



r2{ui,U2;Vl,V2) = J (fzi(fz2T*^"{zi I Ul,W2)£'^t''(^l - ^2)r*'''"'(^2 | Vx,V2) (4.21) 

here the three-line vertices are defined in Eqs. (B.17)-(B.19) and (B.21), D^^,{zi — Z2) = iA|^^,(zi — Z2) with 
Ayt'(^i ~ -^2) defined in Eq. (4.7) and r3(ui, W2; ^'i, ^^2) defined in Eq. (B.22) is the quark-antiquark four-line proper 
vertex. After substituting the expressions in Eqs. (4.18)-(4.21), which are now given in the presence of source J, 
into Eq. (B.6) and completing the differentiation, the one-particle irreducible decomposition of the Green's function 
G'^i^{xi I x\X2;yi,y2) will be found and, thereby, we can write 

2 

Y,Gc'{xi \ xiX2;yi,y2) "^^TG^Axi \xiX2;yi,y2) + ^TGtAx2 \ xiX2;yi,y2) 



^ (4.22) 



3 

= E Gj{xiX2;yi,y2) 



where 



Gi{xiX2;yi,y2) 

= 111 d\d\ E f^r[Au(a^i I xi;«i)5^(x2 - U2) + Sf{xi - ui) (4.23) 
xA^^°-{xi I X2;u2)]T{ui,U2;vi,V2)Sf{vi -yi)S'jr{v2 -1/2) 
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and 



G2{xiX2;yi,y2) 

2 2 

= / n d\d\ E ^TM^i - ui)S^p{x2 - U2)T{ui,U2; ^1,^2) (4.24) 
j=i i=i 

x[A°(a;i I t;i;2/i)S'f,(v2 -^2) + A^«(a;i | V2]y2)SF{vi - yi)] 
Gs{xiX2;y\,y2) 

= / n (i^ujrf^v, E ^7Sf{xi - ui)S^p{x2 - U2) (4.25) 

j=i i=i 

xT''i^{xi I ui,U2;vi,V2)Sf{vi -yi)S^{v2 -1/2) 

where r"''(a;i | ui,U2;vi,V2) is a kind of five-line vertex which is defined in Eq. (B.27) and will be specified soon 
later. 

Before specifying the function r°-'^{xi \ mi, M2; vi, 1^2), we first analyze the expressions in Eqs. (4.23) and (4.24). 
According to the expressions in Eqs. (4.11), (4.13) and (4.18), Eq. (4.23) may be represented as 

Gi{xiX2;yi,y2) = j d'^zid'^Z2[T,{xi,X2;zi,Z2) + Kt{xi,X2;zi,Z2)]Gc{zi,Z2;yi,y2) (4.26) 

where ^{xi,X2; Zi, Z2) and Kt{xi, X2] zi, Z2) wore respectively described in Eqs. (4.12) and (4.14). In view of the 
decompositions in Eqs. (4.5) and (4.8), Eq. (4.24) may be written in the form 

G2{xiX2;yi,y2) = j d^zid^Z2Ki{xi,X2;zi,Z2)SF{zi - yi)Sp{z2 - y2) (4.27) 
where 

2 2 

Ki{xi,X2;zi,Z2) = J II d'^Ujd'^vY^O,^''SF{xi-ui)S^{x2-U2) . . 

j=i i=i (4.^8) 

x[T{ui,U2;v,Z2)T:^{xi \ v; zi) + T{ui,U2; zi,v)J:^''{xi \v;z2)] 

in which T,'^{xi \ v; zi) and T,'^'^{xi \ v; Z2) were represented in Eqs (4.6) and (4.9) respectively. 

Let us turn to the five-line vertex function contained in Eq. (4.25). This vertex is two-particle reducible (or say, 
B-S reducible) although it is one-particle-irreducible. Therefore, it can be decomposed into a B-S irreducible part 
r"^ and a B-S reducible part 

V^^ixi I ui,U2;vi,V2) = T^ji^ixi | ui,U2;vi,V2) + r''jl'p{xi \ ui,U2;vi,V2) (4.29) 

In order to exhibit the above decomposition specifically, as mentioned in Appendix B, we may insert Eqs. (4.19)- 
(4.21), which are now given in the case of presence of the source J, into Eq. (B.27) and complete the differentiation 
with respect to the source J°'^{xi). After completing the differentiations shown in Eqs. (B.23) and (B.25), we can 
write 

3 

T''^'{x^ I ui,U2;vi,V2) = ^r;'^(a^i | uuU2;vi,V2) (4.30) 

where r°^(a;i | ui,U2;vi,V2) are given by the differential of the functions rj(ui, U2; ui, 1^2) in Eq. (4.19) and separately 
shown below 

V°l^{xi I ui,U2;vi,v2) 

= -J d''zd^Zid^Z2Dlll{x, - z){Tl^-\z,Zi I Ui,Vi)Dli,izi - Z2)r^'^'(^2 I U2,V2) (4.31) 
+r-^(^l I Ui,Vi)[Df^y{z, - Z2)Tli-^{z,Z2 I U2,V2) + Iilix'{^.^l.Z2)T-''-\z2 \ U2,V2)]} 

in which besides the gluon-quark and gluon-antiquark three-line vertices mentioned before, there occur the gluon- 
quark four-line proper vertex V^j^^{z, z\ \ u\,v\) and the corresponding gluon-antiquark one r^^_)^(2;, z\ \ U\,v\) which 
are defined respectively in Eqs. (B.28) and (B.29) and 
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nl-^,{z,z^,Z2) = I d^urd^U2D^^,izi - ui)V^^:^:{z,ui,U2)D^J,i{u2 - Z2) (4.32) 
in which r^^,J, {z,ui,U2) is the gluon three- Hne proper vertex defined in Eq. (B.24), 

Tl'^{Xi I Ui,U2\Vx,V2) 

= J d\zd^z,cl^Z2D';i{x, - z){ri-^^*iz,zi I uuU2)Dli,{zi - Z2)r;:'^'*(z2 \VI,V2) (4.33) 

+r'=^*(zi I Ui,U2)[Dli,{z, - Z2)Tlt'^*{z,Z2 I VuV2)+Ul-^\,{z,ZuZ2) T"' ^' * {Z2 \ VUV2)]} 

in which T'^'^'{z, z\ \ wi, U2) and '^^^^iz, z\ \ V\,V2) are the gluon-quark-antiquark four-Hne proper vertices defined in 
Eqs. (B.30) and (B.31) and especially 

Vl^{xi I uuU2\v^,V2) = J d^zDf,{xi - z)T'"'{z I uuU2;vuV2) (4.34) 

in which r'"^(z | ui,U2;wi,f2) stands for the gluon-quark-antiquark five-line proper vertex defined in Eq. (B.33). 

It is pointed out that the vertices formulated in Eqs. (4.31) and (4.33) are not only one-particle-irreducible, but also 
B-S irreducible. This point can be seen more clearly when the vertices are represented by their Feynman diagrams. 
From the diagrams, one can find that it is impossible to divide each of the diagrams into two unconnected parts by 
cutting two fermion lines. However, the five-line proper vertex r'"^(z | ui,U2;vi,V2) in Eq. (4.34) is B-S reducible. 
It can be decomposed into a B-S reducible part T^^^ and a B-S irreducible part F^"^, 

T'"'{z I Ui,U2;Vi,V2) = rj'^j(z I Ui,U2;Vi,V2) +r^£;(z I Ui,U2;Vi,V2) (4.35) 

This enables us to write Eq. (4.34) as follows: 

^Ti^i I Ui,U2;Vi,V2) =r'^l^{Xi I Ui,U2;Vi,V2) +T'^2{^i I Ui,U2;Vi,V2) (4.36) 

where 

Tfiixi I uuU2;vi,V2) = J d^zDf^{xi-z)V%{z \ u^,U2;v^,V2) (4.37) 

and 

Tlt^{xi I uuU2;vuV2) = j d^zDf^{xi - zf^^si^ I uuU2;vuV2) (4.38) 

From the above statement, it is clearly seen that the B-S irreducible part of the vertex in Eq. (4.29) is given by the 
sum 

+Tl^{xi\ui,U2:,Vi,V2) +Tl'({xi\ui,U2:,Vi,V2) ^' ' 



where the three terms on the RHS of Eq. (4.39) were given in Eqs. (4.31), (4.33) and (4.37) respectively. While, the 

B-S reducible part in Eq. (4.29) is given by Eq. (4.38) 

r^'^(Xi I Ul,U2;Wl,W2) = r3^(a;i | Mi,U2;Wi,W2) (4.40) 

Based on the decomposition in Eq. (4.29), the function in Eq. (4.25) will be decomposed into 

G3ixiX2;yi,y2) = G^^\xiX2;yi,y2) + Gf\xiX2;yuy2) (4.41) 
where 

G'i^\xiX2]yi,y2) = j d*zid*Z2K2{xi,X2;zi,Z2)SF{zi - yi)Sp{z2 - y2) (4.42) 
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in which 

/2 2 
Y\_d*Uj'^il.°;'^SF{xl -Ul)Sp{X2 -M2)r^^(xi I Ul,U2;Zl,Z2) (4.43) 
j=l i=l 

and 




(4.44) 



^^'liEi^i I ui,U2;vi,V2)Sf{vi -yi)S^{v2 -2/2) 



It is emphasized that due to the B-S irreducibiUty of the vertex r'^'^{xi \ ui, U2;vi,v2), the function G3 '(x-\_^X2; 2/1,2/2) 
can only be written in the form as shown in Eq. (4.42). While, since the vertex r^'^(xi | ui, U2; wi, ^2) is B-S reducible, 

as was similarly done for the Green's function Q^^^{xi,X2',yi,y2), the function 03^'' (a;i,a;2; 2/1, 2/2), as a part of the 
connected Green's function, must be represented in the B-S reducible form 

G^^\xiX2;yi,y2) = j d'^zid'^Z2K{xi,X2;zi,Z2)Gc{ziZ2;yi,y2) (4.45) 

where K{xi,X2]Zi,Z2) is a kind of kernel needed to be determined later. 



D. Complete expression of the interaction kernel 

Up to the present, the irreducible decompositions of the functions Q'^^^ (2^1, a;2; 2/1, 2/2) (* = 1, 2) appearing in Eq. (2.24) 
have been completed. Now, let us first sum up the expressions given in Eqs. (4.11) and (4.13) which correspond to 

the second and third terms in Eq. (4.4) and the expression in Eq. (4.26) for the fimction G'i(a;i,a;2; 2/1, 2/2) which is 
contained in the connected Green's functions G'^^{xi \ a;i,X2; 2/1, 2/2)- The summation yields 

j d^zid^Z2[T,{xi,X2;zi,Z2) + Kt{xi,X2;zi,Z2)]Q{zi^Z2;yi,y2) (4.46) 

where the relation in Eq. (4.2) has been considered. Then, we combine the expression in Eq. (4.16) which corresponds 
to the last term in Eq. (4.4) and the expressions in Eqs. (4.27) and (4.42) which come from the B-S irreducible part 
of connected Green's functions G"^{xi \ xi, 0:2; 2/1, 2/2) and obtain 



j d'^zid'^Z2[Ks{xi,X2;zi,Z2) + Ki{xi,X2;zi,Z2) + K2{xi,X2; zi, Z2)]Sf{zi - yi)Sp{z2 - 2/2) 



(4.47) 



The final decomposition of the sum of the functions Q'^^^ixi \ xi,X2]yi,y2) and Q^'^^{x2 \ xi,X2',yi,y2) will be given 
by the sum of Eqs. (4.45)-(4.47). Obviously, in order to make the D-S equation to be closed, the kernel In Eq. (4.45) 
must be 

K{xi,X2;zi,Z2) = Ks{xi,X2;zi,Z2) + Ki{xi,X2; zi, Z2) + K2{xi,X2; zi, Z2) (4.48) 
Thus, the summation of Eqs. (4.45) and (4.47) gives 

j d^zid^Z2[Ks{xi,X2;zi,Z2) + Ki{xi,X2;zi,Z2) + K2{xi,X2;zi,Z2)]Q{ziZ2;yi,y2) (4.49) 
Combining Eqs. (4.46) and (4.49), we eventually arrive at 

'^Q^''\xi\ xi,X2\yi,y2) = I d'^zid'^Z2K{xi,X2;zi,Z2)Q{ziZ2;yi,y2) (4.50) 
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where 



K{xi,X2;zi,Z2) = T,{xi,x2;zi,Z2) + Kt{xi,X2;zi,Z2) ,^ 
+Ks{xi,X2;zi,Z2) + Ki{xi,X2;zi,Z2) + K2{xi,X2;zi,Z2) ^' ' 

This just is the kernel appearing in Eq. (2.38). The five terms on the RHS of Eq. (4.51), as respectively shown in Eqs. 
(4.12), (4.14), (4.17), (4.28) and (4.43), are only represented in terms of the quark, antiquark and gluon propagators 
and some kinds of three, four and five-line proper vertices and therefore exhibits a more specific structure of the 
kernel. The equivalence between the both expressions given in the preceding section and this section for the kernel 
will be illustrated in section 6 for the one-gluon exchange kernels. The exact proof of the equivalence has been done 
by the technique of irreducible decomposition of the Green's functions. From the proof, we find that the expression 
described in this section can surely be obtained from the expression given in the preceding section. 



5. Pauli-Schrodinger equation 

As mentioned in Introduction, the D-S equations formulated in the Dirac spinor space may be reduced to an equivalent 
equations represented in the Pauli spinor space with the help of Dirac spinors. Let us start from the equation given 
in Eq. (2.44). The Dirac spinors are defined as [30] 



where are the Pauli matrices, uj = ("p*^ + m^)^/^, U{j^) and V{^) are the positive energy and negative energy 
spinors respectively. They satisfy the orthonornality relations 

u+{-f)u{-f) = v+{-f)v{r) = i 

and the completeness relation 

A+{f)+A-{f) = 1 (5.4) 
where A+("p*) and A~('p*) are respectively the positive and negative energy state projection operators defined by 

A+{^) = U{-f)U+{-f), K-{-f) = V{r)V+{^) (5.5) 
Define 

then, the two fermion spinors can be written as 

Wab(P,r) = Wa{pt)Wb{n) (5.7) 

With this definition, the completeness relation for two fermion spinors can be represented as 

Y,Wab(P,t)W+(P,t) = l (5.8) 

ab 
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Prcmiltiplying Eq. (2.44) with W^f^{P , if ), sandwiching Eq. (5.8) between the kernel K{P, q, k) and the amphtude 
Xp<;{k) on the RHS of Eq. (2.44) and applying the Dirac equation 

h{t)Wa{t) = aco{t)Wa{t) (5.9) 
we obtain 



cd 



(5.10) 



where a, b,c,d= ±1, P = {E, 

CPab{P,q) = W+(P,lf)XpM (5.11) 

Aab{P, t) = E- aujiipt) - boj2{M) (5.12) 

KabcdiP, q, k) = W+(P, t)K{P, q, k)Wcd(P, t) (5.13) 

Eq. (5.10) is a set of coupled equations satisfied by the amplitudes 4>ab{P, q) each of which is represented in the Pauli 
spinor space and of dimension four. In the infinite-dimensional space of the momentum q or k, according to ab = ++ 
and ab ^ H — h, Eq. (5.10) may be, in the matrix form, separately written as 

A++(p)</.++(P) = i^++++(P)</.++(P)+ Yl K++UP)MP) (5.14) 

and 

Aab{P)cl)ab{P) = Kab++{P)<P++{P) + ^ Kabcd{P)ct>cd{P) (5.15) 

cdji++ 

where ab ^ ++ and the terms related to (j)+j^{P) have been separated out from the others. Furthermore, In the 
three-dimensional spinor space spanned by ^ab{P) with ab ^ ++, Eqs. (5.14) and (5.15) may be written in the full 
matrix form 

A+(P)V(P) = K+{P)ij{P) + K\P)ct>{P) (5.16) 

and 

m = G{p)i,{p) + G{p)m (5.17) 

where ^(P) = </.++(P), A+(P) = A++(P), K+{P) = K++++{P), while, </>(P) = {MP)}, k\p) = {X++cd(P)}, 
G{P) = {Kab++{P) / Aab{P)} and G(P) = {Kabcd{P) / Aab{P)} represent the matrices in the three-dimensional 
spinor space. Solving the equation (5.17), we obtain 

4>{P) = x_g(p) g(-P)^(-P) (5-18) 
Substituting the above expression into Eq. (5.16), we finally arrive at 

A+(P)V(P) = F(P)V(P) (5.19) 
where 

V{P) = i^+(P) +7?'(P)^— l^G(P) (5.20) 
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which identifies itself with the interaction Hamiltonian. With the definition 

^ oo 

y3g(p) = EG^"Hn (5.21) 

^ ' n=0 

Eq. (5.20) can be written as 

oc 

1/(P) = ^y(")(P) (5.22) 

n=0 

where 

y(o)(P)=^+(P)_^ 

V^;HP) = AP)G{n_ (5.23) 



Written out explicitly, Eq. (5.19) reads 

The terms in the interaction Hamiltonian in Eq. (5.23) are specified as 

y(°)(P,g,fc) =if++++(P,g,fc) (5.25) 



[E - uiiPt) - t^2{nM{P, Q) = J ^(^' 1, kMP, k) (5.24) 



y«(P,g,fc)= y f dH K++a^iP^q,l)K^,++iPJ,k) 



and so on, where for simplicity of representation, we have defined u)i(T) = u)i{r]i^ + ~t) and W2("^*) = ^2(772^ — 7). 
In the center of mass frame, ZUi( I ) = I ) {i = 1,2). Eq. (5.24) is the equation satisfied by the positive energy 
state amplitude tp{P, q) which is of dimension four in the two-fermion Pauli spinor space. This is the reason why the 
above equation is called Pauli-Schrodinger (P-S) equation. 

By the same procedure, the D-S equation in Eq. (2.45) may also be reduced to a corresponding P-S equation as 
represented in the following 

[90 - V2Mpt) + ViMPimP, Q)= j ^(^' 9' ^)^(^' ^) (5-28) 

where qa is the relative energy and V{P, q, k) is a kind of interaction Hamiltonian which can be written out from the 
expression of V{P, q, k) by the replacement: K{P, q, k) K{P, q, k) and Aab{P, ~'f) Aab(P, q) = qo — a??2t^i(pi) + 
briiui2{p2)- For the three-dimensional D-S equation, the P-S equation in Eq. (5.28) disappears. We are left only with 
a three-dimensional P-S equation derived from Eq. (2.56) such that 

/d^k — » — > 
^V{P,t, k)^{P, k) (5.29) 

where the Hamiltonian V{P,lf,'^) formally has the same expressions as written in Eqs. (5.22) and (5.25)-(5.27) 
except that the four-dimensional kernel K{P, q, k) in those expressions is now replaced by the three-dimensional one 
K{P, "g*, k ) which is the Fourier transform of the kernel in Eq. (3.40) 
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It is worthy to point out that for a given kernel in the D-S equation, there are a series of terms (the ladder 
diagrams) to appear in the interaction Hamiltonian in the P-S equation. If the D-S equation with a given kernel 
could be solved, the contribution arising from a series of ladder diagrams characterized by the series of terms in 
the Hamiltonian are precisely taken into account. Another point we would like to stress is that as seen from Eqs. 
(5.26) and (5.27), the negative energy state only acts as intermediate states to appear in the interaction Hamiltonian. 
Particularly, for the bound state, the positive energy state does not appear in the intermediate states. While, for 
the scattering state P-S equation as discussed in Ref. [30], the intermediate states in the interaction Hamiltonian 
must include the positive energy state. In this case, the series expansion of the interaction Hamiltonian in Eq. (5.22) 
has an one-to-one correspondence with the perturbative expansion of the S-matrix. The above statement reveals an 
essential difference between the interactions taking place in the bound state and the scattering state. 



6. One-gluon exchange kernels 

In this section, wc limit ourself to give a brief derivation and description of the one gluon exchange kernel (OGEK). 
First we discuss the t-channel OGEK and the s-channel OGEK appearing in the four- dimensional D-S equation to 
illustrate the equivalence between the expressions of the interaction kernel derived in the sections 3 and 4. Then, we 
show the OGEK in the three-dimensional D-S equation and the corresponding Hamiltonian in the P-S equation. 



A. The t-channel one-gluon exchange kernel 

The exact form of the t-channel OGEK was represented in Eq. (4.14) with 5]°(a;2 | Xi;zi) and T,'^°-{xi \ X2; z-i) given 
in Eqs. (4.6) and (4.9). In the lowest-order approximation, the propagators and the vertices in Eqs. (4.6) and (4.9) 
are taken respectively to be the free ones and the bare ones. The bare vertices are of the form 

V^^{ux I M2, zi) = -ig-i-TH\ux - W2)<5*(m2 - ^i) 

^\^(u^ I W2, Z2) = -igYT's^ui - U2)6^{u2 - ^2) ^ ' ' 

With the vertices given above, the kernel in Eq. (4.14) becomes 

K^X - Y, X, y) = ig^T-T\^';^,{x2 - yi)SFix, - yihhh25\x2 - 2/2) o) 
+A-t(xi - y2)SUx2 - y2h^7Hnxi - yi)} ^ ■ ' 

Prom now on, the Sf{x — y) and A^^(x — y) in the above are understood to be free propagators. By the Fourier 
transformation, we get in the momentum space that 

K^^{P, q, k) = SiP, q)ig^T'^T' I^fM ' kHl2 (6-3) 
where 

S{P, q) = SF{p,h° + SUP2H = {Sf{pi) + S],{p2)hh2 (6.4) 
in which 

SpM = Sriph" (6.5) 

= ' ^ = . + . (6.6) 

Po-h{p)+is po-u!{p)+is po + u){p) -ts 

here /i(l?) is the free fermion Hamiltonian, A+(j?) and A^('jt) were defined in Eq. (5.5). 

Let us turn to derive the above kernel from the closed expression presented in Eq. (3.30). In the perturbative 
approximation of order g"^, only the first and second terms in Eq. (3.30) can contribute to the OGEK. In the first 
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term which was defined in Eqs. (3.31), (3.23) and (3.21), there are four terms: two represent the self-energies of 
quark and antiquark and the other two are related to the t-channcl one-gluon exchange interaction which is concerned 
here only. The A°(a:2 | xi,zi) and A^"(a;i | X2,Z2) in Eq. (3.21) was respectively represented in Eqs. (4.5), (4.6), 
(4.8) and (4.9). When the vertices are taken to be the bare ones shown in Eq. (6.1), the terms contained in the 
TZ{xi,X2', Zi, Z2) which contributes to the OGEK can be written as 

Hlixi,X2; zi, Z2) = ig^T-T' J dM^f.i^2 - uh^j^-flSpix, - u)j^Sf{u - z,)S^{x2 - ^2) (q 7^ 

+A-t(xi - yH^^j1SUx2 - uh^SUu - z2)6Hxi - z,)} ^ ' ^ 

Next, we turn to the second term in Eq. ( 3.30) which was defined in Eqs. (3.32) and (3.14). Through a perturbative 
calculation ofthe Green's function 5°J^(a;j,2;j | Xi,X2; zi, Z2) defined in Eq. (3.14) or performing a decomposition of the 
Green's function into the connected ones, one may find that there are a function iA'^j^{xi— Zj)SF{xi — zi)Sp{x2— Z2) 

included in the function (?^(^(xi, Zj \ xi, X2] zi, Z2) which just is related to the t-channel OGEK when i ^ j. Thus, 
according to Eq. (3.32), the term included in the Q{xi,X2', zi, Z2) which contributes to the OGEK may be written 
as: 

Hl{xuX2; zu Z2) = ig^T-r'' {Aliixi - Z2hh^SF{xi - z^)S<j,{x2 - ^2)7272° (6 8) 

+AfAx2 - zihh^SUx2 - Z2)SFix^ - ^i)7r7?} ^ " ' 

Substituting Eqs. (6.7) and (6.8) into Eq. (3.30), we will obtain the expression of the t-channel OGEK K^{X—Y, x, y). 
By Fourier transformation, its expression given in the momentum space may be found to be 

2 

K°{P,q,k) = J2HKP,q,k)S-\P,k) (6.9) 



where 



and 

miP, q, k) = tg^T-r'AHiq - fc){7?7r 5^(51)^^(^2)72^72° f6 11) 

Employing the representation of fermion propagator denoted in Eqs. (6.5) and (6.6) and noticing 

S-\P, k) = 7?72°[^F(gi) + SM]-' (6.12) 

one may exactly obtain from Eqs. (6.9)- (6. 11) the expression denoted in Eq. (6.3). Thus, the equivalence between 

the both expressions of the D-S kernel derived in sections 3 and 4 is proved in the lowest order approximation. 

Now let us focus on the three-dimensional t-channel OGEK which was derived for the first time in Ref. [15, 16]. 
For comparison with the four-dimensional kernel, it is necessary to give this kernel a further description based on 
the closed expression formulated in Eqs. (3.40)-(3.47). Analogous to the four- dimensional case, in the lowest order 
approximation, only the first two terms in Eq. (3.40) can contribute to the three-dimensional t-channel OGEK. 
Therefore, we can write 

K°{xi,x2;yi,y2;h -t2) = / d^zid^2;2^-H"i(^,^;^, 22;*! -i2)'S'~^(:zi,:z2;yi,y2) (6.13) 

where Hl{xi, xi; Zi, zi; ti — ^2) arises from Eqs. (3.41) and (3.42) with the three-point Green's functions in Eq. (3.42 
) being given by 

^'^{xi I xi, zi;ti - t2) = -ig J d^uduoAlfl{x2 - - uq) , . 

xSF{xt-lt;h-uo)YT''SFilt -n;uo-t2)} ^' ' 
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A^"(xi I X2,Z2;ti - t2) = -ig J d^uduo^1,i{xi -lf;ti- uq) 

xS'}r{X2 -lt;ti- Uo)YT S'}r{Tt -Z2]Uo- t2)} 



(6.15) 



t2)SF{ 



Xi 



and H2{xi,X2;zi,Z2;ti — ^2) is derived from Eqs. (3.43) and (3.44) when the terms zA°^(xt — ^;ti 

zl;ti — t2)Sp{x2 — Z2;ti — t2) included in Q^'^{xi, Zj \ xt^X2; ^i, 22; E) with i ^ j are taken into account only. By 
the Fourier transformation, it is not difficult to derive the following expression 



K^(P,r,t,E) = J2Ht(P,r,l^,E)S-\'P,t) 



(6.16) 



i=l 



where 



Hl(P,t,t,E) 



(6.17) 



and 



Hl{P,li,li,E) 



(6.16 



+nTA-iik -lf;E-ko-qo)SFm.qo)SU<B,ko)n, } 



The integrals over go and ko can easily be calculated by applying the Cauchy theorem in the complex planes of qo 
and ko- Since QCD is an unitary theory, the matrix element of the kernel between the spinor wave functions is 
independent of the gauge parameter. Therefore, we only need to show the result given in the Feynman gauge. In 
this gauge, noticing the representation of the gluon propagator 



Ql-'^ +ie 

1 



Qo- 



(6.19) 



and the expression shown in Eqs. (6.5) and (6.6), it can be found that 



if - k 

+A-(p5)7072^A+(gJ)] + 



j{P2)+'^{q2)+ 
1 



7?7r[A+(p3)72°72^A-(5J) 



■g*- k 



g*- k 

(5?)]}7?72° 



(6.20) 



and 



if - k 



E-u(p2)-u>(qi)- 



if- k 



+ 



A+(Pl)7i°7r72°72^A+(g2) _ A' (P2 )72°72m7?7;' A' (gl) 



E-w{pi)-L0{q2)- 



if- k 



A-(pT)7i°7r7g72MA-(g2) 



B+a;(Pi)+w(g2) + 



t-k 



E+uj{p2)+ioiqi) + 

}7?72° 



if - k 



(6.21) 



It is seen that H\{P ,lz , E) is actually independent of the energy E. We would like to emphasize that the 
expressions in Eqs. (6.20) and (6.21) can more directly be obtained from Eqs. (6.10) and (6.11) by the following 
integration 
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Hl(P,t,t,E) = J ^'J^Hl{P,q,k),i = 1,2. 



(6.22) 



Now, we discuss the inverse of the function S(P , k ) which is the Fourier transform of the function in Eq. (2.50). 
The equal-time propagators can be defined in such a manner [34] 



Spi'ie' - if) 



-r-T , 



where 



(6.23) 



(6.24) 



With this representation, the function S(P, It) and its inverse, i.e. the three-dimensional counterparts of those in 
Eqs. (6.4) and (6.12) will be written as 



^ > 2i^uj{qt) oj{q^) 



and 



S-\-p,t) = 2ijhl[ 



or^(gi) , M92)i 



t^(9i) ^^iqi) 



(6.25) 



(6.26) 



When Eqs. (6.20), (6.21) and (6.26) are substituted into Eq. (6.16), one may write out explicitly the expression 
of the three-dimensional t-channel OGEK. On inserting this kernel into the first term of the effective interaction 
Hamiltonian denoted in Eq. (5.25) and employing the orthogonality relations of Dirac spinors and the Dirac equation, 

we are led to [15, 16] 



vl''\'P,t,k^) = g^T^TlAit - t;E)U{n)l'(U{qi)Um)l2y.U{qi) 
where U{~^) was represented in Eq. (5.1) and 
A(9^-"r;£) = - 



(6.27) 



"5*- k 



E-u>{pi)-u:{qi)- 



"5*- k 



+- 



(6.28) 



B-w(Pi)-a;(g'2)- 



k 



is just the exact three-dimensional gluon propagator given in the Feynman gauge which is off-shell because the 
energy E is off-shell. It is noted here that the lowest order interaction Hamiltonian 

V^°\'P,t,k^) is only given 

by the fimction i?|(^, ifi k ,E) in Eq. (6.21) because the function H\{P , 'f, k ,E) in Eq. (6.20) gives a vanishing 
contribution to the lowest order Hamiltonian. 



B. The s-channel one-gluon exchange kernel 

The four-dimensional s-channel OGEK was represented in Eq. (4.17). By means of the charge conjugation of the 
quark field, the vertex in Eq. (4.17) can be expressed as 

r'"(^ I ^1,02),. = -C^xT'-iz I ^2,^i)a, (6.29) 

With this relation and the expression shown in Eq. (6.1), the kernel in the lowest-order approximation can be written 
as 
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K°{xi,X2]Zi,Z2)af)pa = ig{A'^l{xi - Zi){nl'^)a^S*p{Xi -0:2)7/3 



(6.30) 



In the momentum space, it reads 

K°{P, q, kUp, = igAf^{P)L'^^{P, qU{C^''T^)p, (6.31) 
where 

L'^^iP, g)a/3 = {^T)c.^S*p{-p2)^fi + (r!r)/3A^>(pi)aA (6.32) 

in which 

S*p{p) = J d^xS*p{x)e-^''=' = Sf{p)C-^ = C-^S^-pf (6.33) 

Now, let us derive the above kernel from the closed expression in Eq. (3.30). Prom the perturbative calculation, 

It can be found that in the lowest order approximation, only the second term in Eq. (3.30) can contribute to the 
s-channcl OGEK because in the perturbative expansion of the Green's function ^^^(a;^, Zj \ xi,X2; zi, Z2), there is a 

term —iA'^'l{xi — Zj)Sp{xi — X2)Sp{zi — Z2) which is merely related to the s-channel OGEK. Thus, the terms in the 
Q.{x\,X2\ z\,Z2) which contribute to the s-channel OGEK, according to Eq. (3.32) can be written as 

2 

H\xx.X2\zx, Z2) = -Y, i^'Jli^i - Zj)nt''S*p{xi - X2)S*p{zi - Z2)Tt^ (6.34) 

Substituting the above expression into Eq. (3.30), in the momentum space, we have 

K°iP,q,k) = H'iP,q,k)S-\P,k) (6.35) 
where 



W{P, q, kUpxs = iAll{P)L''^{P, qU/^L" (P, k)xs (6.36) 
in which L°-^{P, q)ap was given in (6.32) and 

t\P,k)xs = -\Sp{-qi)xr(pTU +%{q2U(pT)rx\ (6-37) 

here 

S*p{q) = J d^xSp{x)e"'=' = C-^S^piq) = Sl{-q)C-^ (6.38) 

In light of the charge conjugation for the 7-matrix and for the propagators shown in Eqs. (6.33) and (6.38), it is 
easy to find 

t^P, kUS-\P, k)x5pa = 9{Cj''T')p, (6.39) 

With this relation, we see, the kernel in Eq. (6.35) is exactly equal to the one written in Eq. (6.31). This gives a 
further proof of the equivalence between the both expressions of the D-S kernel derived in sections 3 and 4. By the 
charge conjugation, it is not difficult to find 

L''^{P,q)ap = gS{P,q)c.pxr{C-^^T"')xr (6.40) 
where 

5(P,g) = Sf(pi)+5^(P2) (6.41) 
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Therefore, the kernel in Eq. (6.31) can be expressed as 

K^iP, q, k) = i5'A«t(P)S(P, qUxACrT'^)xr{CYT^)p. (6.42) 

At the last of this section, we would like to discuss the three-dimensional form of the s-channel OGEK. In accordance 
with Eq. (3.40), this kernel is represented as 



K°{xi,X2;yi,y2;ti - t2) = j(fzid^Z2H%xi,X2;zi,Z2;ti-t2)S ^{zi,Z2; 



yiM) (6.43) 



where iJ*(Si, X2; , I2; — t2) can be written from Eq. (6.34) by setting x\ = X2 = ti and 2;° = 2° = *2 in the 
equal-time frame, that is 



2 



H^xi, xi; zt, zi; h - t2) ^ -^ J2 " ^jlh " i2)f^r^F(^i " ^2)^f(^i " ^2)^/ (6.44) 

In the momentum space, it is of the form 

H'{Pa,~^\E)„pxs = iAli(P,E)L'^f^(P,r)c.pL'''iP,t)xs (6.45) 
where 

L»^(i^, r)a0 = ma^S*p{-t2h0 + {^T)lS'rS*F{tlU (6-46) 

and 

t''(P,t)x5 = -[SU-tl)Xr{^T)r5+S*p{t2)rSiTlT)rx] (6.47) 

which are the three-dimensional form of the functions in Eqs. (6.32) and (6.37). It is emphasized that in Eq. (6.45), 

only the gluon propagator is dependent on energy E, while, the fcrmion propagators arc energy- independent. By the 
same charge conjugation transformations as shown in Eqs. (6.33) and (6.38), one may obtain a kernel similar to Eq. 
(6.42) 

K°{1^, t, E)^pp„ = ig'A';!l(P, E)S{7^, l}Uxr{CrT'')xr{Ci''T'-)p, (6.48) 
which may also be represented as 

i^°(:?, t; E)^0p, = ig^A^iiP, E)S(P, t)a0Xr{rCT'')xr{Cj''T')p. (6.49) 
where 

s(P,r) = s{P,thh2 (6.50) 

which is the three-dimensional form of Eq. (6.4). 

In the P-S equation, similar to Eq. (6.27), the lowest order interaction Hamiltonian given by the kernel in Eq. 
(6.48), according to Eq. (5.25), will be written as 

Vj°\'P,t,t;E) = g^A;l(P,E)U^ipt)+U0in)+{CrT''U{CrT')Mqt)UA^^ (6.51) 

It should be noted that the positive energy state Dirac spinors used here were defined in Eq. (5.1). The negative 
energy state spinor may be given by the charge conjugation relation: V{'^) = CU{jf ) here C = 7^7° [30]. However, 
the matrix C in Eq. (6.51) is defined by C = 17^7'^ [18]. Correspondingly, the charge conjugation relation between 

the spinor wave functions is given by v^{jt) = Cu^{^)'^ where u'^(jf) = u®(p*)+7" with u^ij^) and v^{^) being 
the positive and negative energy spinor wave functions respectively [18] and represented as 

u\-^) = UitW'{t),v'{t) = V{t)x'it) (6.52) 
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here <^*(^) and x*(^) are the spin wave functions and 

Um = \^U{l^),V{lf) = -i/^7V(1?) (6.53) 



m 

Usually, the S-matrix element given by the kernel in Eq. (6.48) is represented by 

T,(P,^,k^) = K-{pt)+uf{n)+K^s(f,t,t;E)^^p,u;^{qt)u^,^{q$) (6.54) 

For later derivation, it is more convenient to use the expression of the kernel written in Eq. (6.49). On inserting Eq. 
(6.49) into Eq. (6.54) and noticing 

< {Pt)+uf {n)+S{^, t)aisxr = -iul' {pt)up (p3) (6.55) 
which is obtained by applying the Dirac equation and 

one can get 

uv,if,k^) = 9^/s.f,{T^,E)u'-m)rT^v-'Kr2)r-mwT\^'m) (6.57) 

This just is the S-matrix element for the one-gluon exchange interaction taking place in the s-channel. It is easy to 
verify that the above matrix element is independent of the gauge parameter. Therefore, we only need to work in the 
Feynman gauge. In this gauge. 

With this propagator, as shown in Ref. [35], by the charge conjugation and Fierz transformation, Eq. (6.57) can be 

represented in the form 

1/2 

T,(Pa,k^)=[^^^^ «yW(:?, g^,fci^)^.,V'r. (6.59) 
where 

(P, t,k^) = ^'^f' U{pt)U{m^i2Uiqt)U{q^) (6.60) 

£"^-7 +ie 

is the interaction Hamiltonian occurring in the P-S equation in which the spinor is still defined in Eq. (5.1), C's is 
the color matrix 

Cs = ^(A? - AD' (6.61) 
with Xf being the Gell-Mann matrices, Fg is the flavor matrix which has an expression for flavor SU(2) such that 

Fs = ^{1 - Tt ■ T^) (6.62) 
here ft are isospin Pauli matrices and 

ri2 = -/1/2 + 7hl - ^7r72M + ^(7hr)(7l72^) (6.63) 

In the end, it is pointed out that since the matrix clement of the color operator Cs in the qq color singlet vanishes, 
the s-channel OGEK contributes nothing to the qq bound states. However, for many-quark-antiquark systems such 
as TTTT, KK, ttN, KN systems and etc., the contribution of the s-channel OGEK is not negligible and plays an 
important role to the interations taking place in those systems. [36- 38]. 
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7. Discussions and remarks 



In this paper, the D-S equation satisfied by the qq bomid states has been derived from QCD and, especially, the 
interaction kernel in the equation has been given two equivalent closed expressions which were respectively derived 
by making use of the equations of motion obeyed by the Green's functions and the irreducible decomposition of the 
Green's functions. Since the B-S equation is commonly viewed as the correct equation for the bound state problem, 
it is natural to ask what is the relation between the D-S equation and the B-S equation? As shown in Ref. [5], the 
B-S equation may be derived from the D-S equation. In fact, when applying the Lehmann representation in Eq. 
(2.36) to Eqs. (2.16) and (2.17), by the same procedure stated in section.2, one may obtain two D-S equations as 
follows 

{idxi-mx)xp<;{xi,X2) = l d^yid'^y2Ki{xi,X2\yi,y2)xp<;{yi,y2) (7-1) 



{idx^-m2)xp<i{xi,X2) = j d'^yid'^y2K2{xi,X2;yi,y2)xPi{yi,y2) (7.2) 
where 

Ki{x-i,X2;yi,y2) = 7°ii'i(xi,a;2;yi,y2), i = 1,2 (7.3) 

in which Ki{xi,X2;yi,y2) was given in Eq. (3.28). The D-S equations shown in Eqs. (2.38) and (2.39) may directly 
be written out from Eqs. (7.1) and (7.2). 

Operating on Eq. (7.1) with {18x2—^2) or on Eq. (7.2) with {idxi—mi), we have 

{idxi-mi){idx2-m2)xp^ixi,X2) = j d^yid^y2KB{xi,X2;yi,y2)xp<;{yi,y2) (7.4) 

where 

KB{xi,X2;yi,y2) = {idxi-m.i)K2{xi,X2;yi,y2) = {idx2-m2)Ki{xi,X2;yi,y2) (7.5) 

is the B-S interaction kernel whose explicit expression was derived in Ref. [5]. Acting on Eq. (7.4) with the inverse 
of the operator {idxi—mi){idx2—'iTi2), the B-S equation will be recast in an integral equation 

XPi{xi,X2) = j d'^yid'^y2d'^zid'^Z2S]p\xi - zi)S'^°\x2 - Z2)KB{zi,Z2;yi,y2)xP';{yi,y2) (7.6) 

where Sp '' {xi — Zi) and Sp°^ {x2 — Z2) are the free propagators of quark and antiquark respectively. In the momentum 
space, it becomes 

XpM = S^^\p^)Sf\p2) I -0^KB{P,q,k)xP,{k) (7.7) 

Conversely, if we act on Eq. (7.6) with the operators {18x2—^2) and (idxi—mi) respectively, the D-S equations 
in Eqs. (7.1) and (7.2) will be recovered. This shows that to get the B-S equation, we need merely to consider 
the D-S equation. It should be noted that the four- dimensional B-S equation can not directly be transformed into 
the three-dimensional D-S equation. In order to obtain a three-dimensional equation from the four-dimensional B-S 

equation, it is necessary to introduce a certain constraint condition on the relative time (or relative energy) as was 
done in an approximate manner such as the instantaneous approximation or the quasipotential approaches [7-12]. 

As stated above, the four- dimensional D-S equation and the corresponding B-S equation can be derived from one 
another. But, this does not mean that the D-S equation and the B-S equation are fully equivalent to each other, 
similar to the Dirac equation and the K-G equation which can also be derived from each other. As seen from Eq. 
(7.4), the B-S equation is a kind of second-order differential equation in the position space. Therefore, a solution 
to the equation depends on not only the amplitude at time origin, but also the time-differential of the amplitude at 



30 



the time origin as in the case' for K-G equation. Tliis probably is the origin that causes the B-S equation to have 
the unphysical solutions of negative norm. In order to exclude the unphysical solutions, as mentioned before, the 
common procedure is to recast the four- dimensional B-S equation in a three-dimensional form by eliminating the 
relative time (or the relative energy) from the equation. Certainly, the three-dimensional equation, particularly, the 
exact version of the equation presented in the sections 2 and 3 is much convenient to use in solving the eigenvalue 
problem. However, since the three-dimensional equation loses the Lorentz-covariance of a relativistic dynamics, it 
is sometimes not suitable for carrying out extensive theoretical analyses, for instance, to perform the irreducible 
decomposition of the Green's functions contained in the kernel given in Eq. (3.30). The decomposition can readily 
be done in the four-dimensional form as shown in section 4. At this point, we may ask whether the relativistic 
bound state problem can be solved Lorentz-covariantly in the Minkowski space without occurrence of the unphysical 
solutions? The answer should be positive because the Lorentz-covariance of the equation implies that one may work 
in any Lorentz frame and gets the same result. Let us turn to the D-S equations shown in Eqs. (2.41) and (2.42) 
which are represented in the position space. Clearly, the equation in Eq.(2.42) is a first-order differential equation 
of Schrodinger-type. One may first solve this equation to get an amplitude which describes the evolution of the 
amplitude with the relative time t. and then substitute this amplitude into Eq.(2.41) to solve the eigenvalue E and 
the amplitude Xp?(^)- In solving these equations, we only need the initial conditions of the amplitude at the time 
origin without concerning the initial conditions of the time-differentials of the amplitude. Therefore, the unphysical 
solutions would not appear in this case. In this sense, we can say, the D-S equation derived in this paper, as it provides 
a new formulation of the relativistic equation for the two fermion bound system, gives a suitable prescription to solve 
the four-dimensional equation. Moreover, based on the relation denoted in Eq. (7.5), the B-S kernel may conveniently 
be evaluated from the D-S kernel. In comparison with the closed expression of the B-S kernel derived in Ref. [5], the 
D-S kernel shown in Eq. (3.30) is rather simpler. The main contribution to the D-S kernel is given by the Green's 
function Gfiti^i^Uj I Xi,X2;yi,y2) written in Eq. (3.14). While, the B-S kernel concerns more complicated Green's 
functions such as [5] 

0fj\Ti^i'^2,yi,y2 I xi,x2;yi,y2) 



= 0+ 



(7.8) 



T{N[A-{x^)Ai{x2)Hxi)rix2)]N[Al{y,)Ai{y2my,)^ (t/2)] 

which gives the major contribution to the B-S kernel. In particular, in comparison of the four-dimensional kernel 
represented in Eqs. (3.30)-(3.33) with the three-dimensional counterpart written in Eqs. (3.40-)-(3.43), we see, there 
is an one-to-one correspondence between the both kernels. Therefore, to calculate the three-dimensional kernel, 
one may first c;aleulate the four-dimensional one and then convert it to the three-dimensional form according the 
correspondence relation between the both of them. Since the four-dimensional D-S equation is Lorentz-covariant, its 
kernel can conveniently be analyzed and calculated by means of the familiar technique developed in the covariant 
quantum field theory. 

In the end, we would like to address that unlike the Dyson-Schwinger equation [39, 40] which contains an infinite 
set of equations, the D-S equation derived in this paper is of a closed form with a closed expression of the kernel as 
given in section 3 or section 4. The kernel can easily be calculated by the perturbation method. For example, in 
the perturbative calculation of the kernel given in section 3 which contains only a few types of Green's functions, 
we only need the familiar perturbative expansions of the Green's functions without concerning the calculation of 
other more-point Green's functions as it is necessary to be done for the Dyson-Schwinger equation. Especially, each 
of the Green's functions can be represented in the form of functional integral and is possible to be evaluated by a 
nonperturbative method as suggested, for example, by the lattice gauge theory. Therefore, the expression of the 
kernel given in this paper provides a new formalism for exploring the QCD nonperturbative effect and the quark 
confinement which are important for the formation of a qq bound state. In the ordinary quark potential model [3, 41], 
the quark confinement is usually simulated by a linear potential which was suggested by the lattice computation of a 
Wilson loop and by the area law [42-44] . Obviously, this simulation is oversimplified. For the purpose of investigating 
the quark confinement, it is appropriate to start from the kernel given in this paper for the case that the quark and 
the antiquark have different flavors. In this case, all the Green's functions become the ordinary ones as shown in Eqs. 
(2.5), (2.13) and (3.11). Since the kernel derived in this paper contains all the interactions taking place in the bound 
state and includes the color-spin matrices $1°'' defined in Eq. (2.23) in it, it is anticipated that a nonperturbative 
calculation of this kernel would give a sophisticated confining potential which includes not only its spatial form, but 
also its spin and color structures. This just is the advantage of the formalism of D-S equation presented in this paper. 
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9. Appendix A: Derivation of equations of motion satisfied by the Green's functions 

This appendix is used to derive the equations of motion satisfied by the quark-antiquark two and four-point Green's 
functions. These equations may be derived from the following QCD generating functional [18]. 

Z[J, % V,l^] = j^j T^iA V^, V, C, C)e'' (A.l) 
where 

I = J d^x[C + r^'Al + ri^+^r] + lC + C^] (A.2) 
in which £ is the effective Lagrangian of QCD 

£ = li^{id-m + gA)i, - ^F'^^'^F^;, - ^(^^^^^ + C''d'^{D;'>C>') (A.3) 
here A ='y^T'^A^ with being the vector potentials of gluon fields, 

F;, = d^A- - d.A- + gr'"'AlAl (A.4) 
are the strength tensors of the gluon field, 

^ab ^ gabg^ _ gpbc^c^ (A.5) 

are the covariant derivatives, TT ^C^ represent the ghost fields, and J"^'^ ,ri,T],^ and denote the external sources 
coupled to the gluon, quark and ghost fields respectively. By the charge conjugation transformations shown in Eq. 
(2.2) for the quark fields and in the following for the external sources 

rj" = Crf.rf = -ri^C'^ (A.6) 
it is easy to prove the relation 

ip{id-m + gA)ip + TyV + V'»? = i^^iid-m + gA^-il^" + rj^ip" + ij^r]" (A.7) 
where A = 7''T°A» . 

A. Equations of motion with respect to the coordinate xi 

Upon taking the functional derivative of the generating functional in Eq. (A.l) with respect to the field function 
V'c(a;i) and considering 

(A.8) 



and 

SI 



^-r , = Vaixi) + [{id^^ - mi)aj + gA{xi)aj]'ipj{xi) (A.9) 
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it can be found that 

where the fields A'^^{xi) and tl)~f{xi) have been replaced by the derivatives of the generating functional with respect 
to the sources J°'^{xi) and Ti^{xi) and each of the subscripts a, (3 and 7 marks the components of color, flavor and 
spinor. Differentiating Eq. (A. 10) with respect to the source r}p{y{) and then setting all the sources to vanish, we 
obtain the equation satisfied by the quark propagator [18] 

[{idxi - mi + S)S'ir]„p(a;i, yi) = 5ctp5'^{xi - yi) (A.H) 

where 

{T,SF)ap{xi,yi) = J d'^ziT,{xi,zi)ajSF{zi - yi)jp . 

= {T-n^^A-ix,\xuyi)-y0 ^'^■'^> 

here T,{xi,zi) stands for the quark proper self-energy and A."{xi \ xi,yi)jp was defined in the first equality of Eq. 
(3.11). 

Let us turn to derive the equation of motion satisfied by the Green's function defined in Eq. (2.1). In doing this, 
we need first to derive the equations of motion obeyed by the Green's function defined in Eq. (2.5). By successively 
differentiating Eq. (A. 10) with respect to the sources ^^(2^2), ?7p(yi) and r]'^{y2), noticing the equality in Eq. (A. 7) 
and the following nonvanishing derivatives 

Igg = C.,SHx -y),^ = iC-^US\x - y), 

S0i = iC-%,SHx-y)/^ = C^,5\x-y) ^ ' > 

we have 

+'?"(-^i)7W^TO^7rw^ " [(^^-^ - + (r"'^)«7«7^] (A.14) 

When all the sources are set to be zero, one immediately derives Eq. (2.10) from Eq. (A.14). 

It is noted that the first equation in Eq. (2.14) may directly derived from Eq. (A. 11) by the charge conjugation 
transformation represented in Eq. (2.2) or by differentiating Eq. (A. 10) with respect to the source rf^{x2)- 

B. Equations of motion with respect to the coordinate X2 

When taking the derivative of the generating functional in Eq. (A.l) with respect to the field variable ip''^{x2) and 
noticing the relation in Eq. (A. 7), by the same procedure as described in Eqs. (A.8)-(A.10), one may obtain 

W,ix2) + [i^^.2 - m2h. + = (A.15) 

Differentiating Eq. (A.15) with respect to r}%{y2) and then setting all the sources to be zero, one can get the equation 
for the antiquark propagator 

[{id^2 -m2 + S^)5^];3,(ar2, y2) = 50„5\x2 - ^2) (A.16) 
where 



{T.^S^p)p,{x2,y2) = J d^Z2^^{X2,Z2)0xS''p{Z2-y2)xa 
= {T'"')/3xA^J'{x2 I X2,y2)xa 
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here E'^(x2,22) is the antiquark proper self-energy and kf{x2 I 2^2, 2/2) Act was represented in the second equality of 
Eq. (3.12). 

Similarly, Upon differentiating Eq. (A. 15) with respect to the sources ri^{xi), rip{y\) and 77^(2/2), one gets 



+'0''^ hs^^ix,)sSyi)Sv^AV2) ~ [(^^^^ " '^2)/3A + (r''")/3A ,^jb^(^^) ] (A.18) 

^ \Z = 

Sri^{xi)S7]l{x2)Srip(yi)Sr]^(y2) > 

Setting all the sources to vanish, we directly obtain Eq. (2.11) from the above equation. 

It is mentioned that the second equation in Eq. (2.14) may directly be derived from Eq. (A. 16) by the charge 
conjugation transformation or by differentiating Eq. (A. 15) with respect to the source r]^{xi). 

C. Equations of motion with respect to the coordinate yi 

By taking the derivative of the generating functional in Eq. (A.l) with respect to the field variable tjjpiyi), following 
the same procedure as deriving Eq. (A. 10), one may get 

On differentiating the above equation with respect to r]^{xi) and then turning off all the sources, we arrive at 

[Spiidy^ +mi -T,)]cp{xi,yi) = -6apS'^{xi - yi) (A.20) 
where 

{SF'^)ap{xi,yi) = Jd^ZiSpixi - Zi)aT'^{zi,yi)rp l\ 0^\ 

= A-{y,\x^,y^U{T'^nrp ^^'^^^ 

If we differentiate Eq. (A. 19) with respect to 7/^(2/2), after letting the sources to be vanishing, we get an equation 

satisfied by the propagator Sp{yi — r/2) as written in the first equation in Eq. (3.3). 

Now let us differentiate Eq. (A. 19) with respect to ^„(a;i), ^^(2:2) and 77^(2/2) and then set all the sources but the 
source J to be zero. By these operations, we get 

G(a;i,a;2; yi,y2)if}rai0 vi + ^^ijrp = -^ap^'^ixx - yi)S^p{x2 - y2)p„ 

-CpX{yi-y2)SUxi-x2)ip + Gl{yi\xuX2\yi,y2)UrA^''nrp ^ ■ ' 



where 



(^[XiX2;yi,y2) = 7=7 — tt-^t — tt-t — rr— 7 — r 77=r='7='7"=o (A.23) 
Sr]{xi)d'n''{x2)dr]{yi)6r]''{y2) 



iS%{x2)Sr]^{y2) 



S* ix. - x.V — ^^^t-^'^'^'^'^] I - IK 25) 



and 



G^iyi I 0:1x2; 2/1,2/2)^ = .^j„^^^^ G(a;i,a:2;2/i, ^2)'^, i = 1,2 (A.26) 
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Once we set J = 0, Eq. (A. 22) will give rise to Eq. (3.1). Furthermore, if we differentiate Eq. (A. 22) with respect 
to J°'^{xi) and subsequently set J = 0, noticing 

GfAxuVj I xiX2;yi,y2) = ^JJ^^r(^;)^^j6^^(^i.^2;yi,y2)'' |j=o (A.27) 

the equations in Eq. (3.8) will immediately be derived. 



D. Equations of motion with respect to the coordinate 2/2 

To derive the equations of motion with respect to 2/2, we need to differentiate the generating functional with respect 
to the field ^^(2/2)- By the same procedure as formulated in Eqs. (A.8)-(A.10), we get 



{^M + -lz^Midy,+m2)sa - (r '^)5.777wTt]}^ = « (A-28) 



Then, the differentiation of the above equation with respect to r]'^{x2) with setting all the sources to vanish subse- 
quently will lead us to 

[S'p^ity, +1712- E^)]/3,(X2, 2/2) = -S(}„S\X2 - y2) (A.29) 

where 

{Sp'2,'')i3a{x2,y2) ^ j d'^Z2S''p{X2- Z2)l35^''{Z2,y2)5a . 

= ^f{y2\x2,y2hs{t'')e. ^ 

Upon differentiating Eq. (A. 28) with respect to r]p{yi) and then turning off all the sources, one can obtain the second 
equation in Eq. (3.3) for the propagator Sp{yi — 2/2) 

Let us differentiate Eq. (A. 28) with respect to r]^{xi),r]1^{x2) and rip{y\) and set all the sources except for the J 
to vanish. As a result, we get 



G{xiX2;yi,y2)i0ps{idy2 +m2)sa- = -S0aS^{x2 - y2)SF{xi -yi 



-CpaS'^iyi - y2)S*p{xi - X2)ip + G^(y2 | xi,X2]yi,y2)ippsi^^'')5a 
where 



'ap 



(A.31) 



and the other Green's functions in the presence of source J were defined before. When we set J = 0, Eq. (A.31) 
straightforwardly yields Eq. (3.2). Finally, on differentiating Eq. (A.31) with respect to J°'^{xi) and subsequently 
setting J = 0, Eq. (3.9) will be derived. 



10. Appendix B: One-particle irreducible decompositions of the connected Green's functions 

Let us begin with the relation between the generating functional for full Green's functions Z[J,r],ri,^,^] and the one 

for connected Green's functions W[J,rj,ri,£^,£] [18, 33] 

Z[J,ri,v,l^] =eMiW[J,V,V,lm (B-l) 

Taking the derivatives of Eq. (B.l) with respect to the sources rj{xi),rj'^{x2),rj(yi) and ?7^(y2) and then setting all 
the sources except for the source J to be zero, one may obtain the following decomposition 
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G{xiX2;yi,y2Y_= Gc{xiX2]y 1,1)2)'' + Sf{xi - yiys^{x2 - 2/2)'^ 
-S*p{xi - X2ys*p{yi - j/2) 



where G{xiX2-.yi,y2y', Sp{x2 - y2)'' , Sp{xi ~X2Y and Sf{xi -yiY were defined in Eqs. (A.23)-(A.25) and (A.32) 
respectively, wliile, Gc{xxX2\yi,y2y and Sp{yi — 1/2)'' are defined by 

Gc{xlX2;yi,y2y = ^ s ^c . \ ^=5)'==r;=r,= =0, (B.3) 

dr]{xi)6r]''{x2)or]{yi)Sri''{y2) 



and 



PSri{yi)Si]''iy2) 



When we set J = 0, Eq. (B.2) will go over to the decomposition shown in Eq. (4.1). Differentiating Eq. (B.2) with 
respect to the source J°'^{xi), we have 

G1,{xi I xiX2;yi,y2y = G'^^{xi \ xiX2;yi,y2y 

+A«(a;i I xi;yiyS^{xi - yiY + Sf{xi - ViY \ 0:2; 2/2)-^ (B.5) 
-h.l*{xi \xi,X2yS*p{yi-y2y-S*p{xi-X2yJS^''*{xi | 2/1,^2)-^ 

where G1^{xi \ a;i, 0:2; 2/1, 2/2)"^ was defined in Eq. (A. 26) with 2/i being replaced by Xi and the other functions are 
defined by 

Gciii^i I xi,X2;yi,y2y = .^j^^^^^ Gc{xiX2;yi,y2y (B.6) 

^^^^""^ I = ilT^^^^^^ - ^""-^^ 

Al^ix, I x,;y,y = ^^^S^x, - y,y (B.8) 



iSJ^f'ixi) 
5 



A-(x. I XI, ^2)^ = -FT^-T^M^i - x^y (B.9) 



and 

Upon setting J = 0, Eq. (B.5) immediately gives rise to the decomposition in Eq. (4.3). 

Now, let us proceed to carry out one-particle-irreducible decompositions of the connected Green's functions on the 

RHS of Eq. (4.4). The decompositions are easily performed with the help of the Legendre transformation which is 
described by the relation between the generating functional of proper vertices F and the one for connected Green's 
functions W [18, 33] 

T[A1,^,^l^,C^,G-\=W[J,%in,l,^\- J d''x[r''A; + r}^ + ^r, + ^C + C^] (B.ll) 
and the relations between the field functions and the external sources 

Mx) = ——,'0 X = -^-r^,Al{x) = — — —^Clx) = — ,C (a;) = - , , (B.12) 
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r]{x) = -JJ—,r]{x) = ■j^,J?,{x) = -TT^,r(x-) = - -a = 77?^ (B.13) 

where the field functions in Eq. (B.12) are all functionals of the external sources in Eq. (B.13) and, simultaneously, 
the sources in Eq. (B.13) are all functionals of the field functions in Eq. (B.12). 

Taking the derivative of the both sides of the first equality in Eq. (B.12) with respect to '^(j/) and employing the 
first relation in Eq. (B.13), one may get 

a z = — ,,,,,, — d, z——^———^ = —0 [x — y) B.i4) 

6^{y)6i){z) Sr]{z)Srj{x) J 8r}{x)5r]{z) 5^{z)5^{y) ^ ' ^ ' 

where we only keep the term on the RHS of Eq. (B.14) which is nonvanishing when the sources arc set to vanish. In 
order to find the one-particle-irreducible decomposition for the quark-gluon three-point Green's functions, one may 
differentiate Eq. (B.14) with respect to the source J°-'^{xi) and then using Eq. (B.14) once again. By this procedure, 
it can be derived that 

- J a ^« «1« "2^j„^(^.)5j6.(„^) STj{xi)Sr,(u2) /p. T 

6^w [D.io) 



where the coordinates in Eq. (B.14) have been appropriately changed. When all the sources are set to be zero, 
noticing the definitions given in Eq. (A. 32) where the Z is replaced by iW and in Eq. (B.7) as well as 



r'"^(ni I U2, z) = i U=v^=^=o (B-17) 

the decomposition shown in Eqs. (4.5) and (4.6) straightforwardly follows from Eq. (B.15). Analogously, if we 
replace ri{xj) and ri{yk) by r]'^{xj) and r)''{yk) in Eq. (B.15) and noticing 

the decomposition shown in Eq. (4.8) and (4.9) will be derived. This decomposition may also be derived from Eq. 
(B.15) by the charge conjugation transformation for the quark fields. By this transformation, one may readily derive 
from Eq. (B.15) the decomposition denoted in Eq. (4.15) in which the gluon-quark-antiquark vertex is defined by 

^ ^'^''^''^=' 6At{z)6i.{z,)5nz2) ^""-''^ 

The one-particle-irreducible decomposition of the connected Green's fimction Gc(a;i X2; j/i, 1/2) can be derived by 
the same procedure as obtaining Eq. (B.15). On differentiating Eq. (B.14) with respect to r]^{x2) and rf{y2) and 
setting all the sources but the source J to vanish, one may obtain 

2 

Gc{xiX2;yi,y2y = J JJ d'^Uid'^ViSpixi - uiY S'}r{x2 - U2y (B 20) 

xr(ui,u2;t^i,t^2)'^<S'F(wi -yiys%{v2 -y2y 

where the four-point connected Green's function and the propagators given in the presence of the sources were 
defined before and the function r(Mi, U2; wi, W2)"^ is formally the same as that defined in Eqs. (4.19)-(4.21). When 
the source J is turned off, Eq. (B.20) directly goes over to the decomposition in Eq. (4.18) with the vertices in Eqs. 
(4.19)-(4.21) being defined in Eqs. (B.17)-(B.19) and in the following 
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X3-r 

which is the charge conjugate to the vertex T^"*{z \ zi,Z2) as well as 



6tl){ui)5lj?'{u2)Stp{vi)6il)''{v2) 



T^iui, U2; vi, V2) = i ^^^ ^^^^ ^ l^^^^-^c^^c^^ (B.22) 



which is the quark-antiquark four-line proper vertex. It is emphasized here that the decomposition of the function 
Gc{xi.X2', jji, 112) in the absence of the source J has the same form as that given in the presence of J . This is because 
the Green's function is defined only by the differentials with respect to the fermion fields as indicated in Eq. (B.3). 

The one-particle irreducible decomposition of the Green's function G'^^{xi \ x\X2; 2/1,2/2) may be derived by starting 
from the expression given in Eq. (B.15) with j, k = \. By differentiating the both sides of Eq. (B.15) with respect to 
the sources r]'^{x2) and 77*^(2/2) and then turning off all the external sources, one may obtain the decomposition of the 
function G'^^{xi \ X1X2] 2/1,2/2) as shown in Eqs. (4.22)-(4.25). Alternatively, the decomposition may also be obtained 
by starting with the expression written in Eq. (B.20). Substituting Eq. (B.20) into Eq. (B.6), then completing 
the differentiation with respect to the source J"'^{xi) and finally setting the source to vanish, one may also derive 
the irreducible decomposition of the function G'^^{xi \ .ti,.T2; j/i, 2/2)- In doing this, it is necessary to perform the 
differentiations of the fermion propagators with respect to the source J'^'^{xi) as shown in Eqs. ( B.7)-(B.10) and 
use their decompositions presented in Eqs. (4.5)-(4.9). In addition, we need to carry out the differentiations of the 
gluon propagator and some vertices with respect to the source J°'^{xi) as shown below. For the gluon propagator 
defined in Eq. (B.16), from its representation in presence of the external source J ,in the same way as deriving the 
decomposition represented in Eqs. (B.15), ( 4.5) and (4.6), one may obtain the one-particle irreducible decomposition 
of the gluon three-point Green's function as follows: 

Kl%{xi,z^,Z2) = —J---A^fMi - Z2y |j=o= / d^zD'^Jlix. - z)n%'-''iz,z,,Z2) (B.23) 



(5J«^(a;i)) 

where Z?«^,(,x, - z) = iA»l{xi - z) and n^^/'"(-2, ^i, Z2) was represented in Eq. (4.32) with the gluon three-line proper 
vertex defined by 



' 5Al{z)5A^i)5Ai{u2) 



TZ2{z,u,,U2) = , U=o (B.24) 



For a proper vertex r„(2i, 2:2, • • •) with a marking the other indices, its derivative with respect to the source J"'^(xj) 
can be represented as 

ra(^i,^2,- • y |7=o= J d^zD^lixi - z)ri^{z,zuZ2,-- (B.25) 



iSJ'^i^{xi)) 
where 

ri^{z,z,,Z2,---) = ^:4|(^I'«(^i' ^2, •••)■' |j=o (B.26) 

According to the procedure stated above, it is not difficult to derive the expressions described in Eqs. (4.30)-(4.33). 
In the expressions, the vertices are defined as follows: 

r ^{x^ \ ui,U2;vuV2) = ^SJ^^^ 1^=0 (^-27) 

T'll(z,zi\ui,vi) =i = L T / n (B.28) 

' ^ SA'"^{z)6A-^{zi)d^iui)Si;{vi) '^=V'=V'=o ^ ^ 
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SAb-{z)dA<=^{zi)Silj\ui)d^^{vi) '^-^ 

—5c* (5^r 

I U,,V,) = *5^,.(,)^^eA(,^)5^(^^)5^c(^^) |a=^=^»=0 (B.31) 

Particularly, by the following differentiation 

it is easy to give the expression in Eq. (4.34) in which 



V' {u 

is the gluon-quark-antiquark five-line proper vertex 
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